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ABSTRACT 


When first introduced to orbital mechanics, students often experience 
difficulty in visualizing a satellite's actual path through space. The Satellite 
Maneuver Evaluation Tool (SMET) seeks to alleviate that problem, SMET is a 
three-dimensional color graphics simulation of satellites in flight. It allows a 
student to interactively modify a satellite's orbital parameters and see the effects 
as the satellites' positions are updated continuously with respect to time. The 
user can change parameters by defining a maneuver or by directly entering a 
change through keyboard or dial inputs. SMET offers the user the capability to 
demonstrate difficult concepts, as well as a method to simulate actual satellite 
maneuvers. Instructors can videotape SMET sessions for classroom 
demonstrations. Instantaneous images can also be saved for redisplay or for 
printout. 

This thesis provides a background on the mathematical formulas modeled by 
SMET. It also incU.des a detailed user's guide for SMET. 
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I. INTRODUCTION 


Students in the Space Systems Operations and Space Systems Engineering 
curricula at the Naval Postgraduate School are required to develop a thorough 
understanding of orbital mechanics. They use equations, tables of data and two- 
dimensional static illustrations to see the relationship between a satellite's orbital 
parameters and its path in space. These students, most of whom arrive with little 
or no background in the field, often have difficulty in visualizing the actual path 
of a satellite’s flight. There is no readily available tool assist in learning orbital 
mechanics; a tool which would allow a student to interactively modify a satellite’s 
orbital parameters and see the effects. 

A mathematical relationship can be represented as two-dimensional graph. A 
graph conveys information more readily to human beings than does a table of the 
underlying numbers. But if more than two variables can be manipulated such 
representations quickly loose there utility. Adding a third dimension to a graph 
and varying it with time can be used to portray more information. But that 
format may not be intuitive to the user. Visual simulation is another approach to 
be considered. Visual simulation is the creation, by computer, of a realistic, 
easily modified, moving image from the mathematical model of a phenomenon. 
Realism implies high-resolution, color graphics. [Ref. l:p. 3-5] 

The Satellite Maneuver Evaluation Tool (SMET) is a three-dimensional color 
graphics simulation of satellites in flight. It offers the user three ways of 
observing a satellite’s position with respect to the Earth. In one window, the 
user is placed in space at an arbitrary point of his choice. From this position, he 
sees the motion of the satellite around a rotating model of the Earth. Another 
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window allows the user to see that portion of the Earth which is visible to the 
satellite. A third window shows a Mercator projection map of the Earth with a 
plot of the satellite's ground track. After a user enters a satellite's parameters, 
the satellite's position is continuously updated with respect to time in all three 
windows. 

Although SMET allows the user to change parameters through the use of 
dials or keyboard input, its greatest utility is in its ability to let the user execute a 
satellite maneuver at any given moment. A maneuver is entered as a velocity 
change (magnitude and direction) or as a combination of three maneuver 
components (flight path angle change, heading angle change and velocity 
magnitude change). The user can also enter the burnout parameters of a launch 
vehicle in order to place a satellite into its orbit. 

The user is also able to read a satellite's position and velocity continuously in 
different coordinate systems and unit types. By allowing the user to enter as 
many satellites as he desires, SMET also lets him visualize the relationship 
between several satellites. As an example, the user could enter several identical 
satellites but vary one of the parameters to visually compare the effect of that 
variation. SMET also allows students to create illustrations and to program 
simulations. 

SMET is written for any Silicon Graphics IRIS 4D/GT, GTX or VGX 
workstation. A personal computer implementation was considered and rejected 
for two reasons. One, there is no commonly available set of three-dimensional 
color graphics routines for personal computers that offers the functionality to 
generate the images desired for SMET. Second, graphics workstations are 
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becoming more affordable all the time. An IRIS is not an extravagant piece of 
hardware. Currently, there are at least 10 IRIS workstations around the campus. 

The next three chapters provide a background on the mathematical equations 
used in SMET’s simulation of satellite motion, maneuver and launch. The 
equations are not derived, in most cases, because that was considered outside of 
the scope of this document. Those derivations can be found in the refe’'''nces or 
other publications listed in the bibliography. All angles in trigonometric 
functions are in radians. Inverse trigonometric functions may require additional 
quadrant checking that is not shown. Although distances and velocities can be 
displayed in metric or english units, SMET actually stores values and performs 
all computations in canonical units. This, however, is transparent to the user. 
Chapter V provides detailed instructions on how to use SMET. Chapter VI 
discusses ways to modify the source code for specific simulations not currently 
modeled in SMET. Appendix A is a consolidated list of common variables in 
orbital mechanics. Appendix B has procedures for recording SMET sessions on 
video tape. It also provides directions on how to print out Postscript files 
generated by SMET. Appendix C is the source code listing for SMET. 





II. SATELLITE MOTION 


A. CLASSICAL ORBITAL PARAMETERS 

A satellite’s path in space can be described by five independent quantities. 
These five quantities, called "orbital elements" or parameters, are sufficient to 
completely describe the size, shape and orientation of the orbital ellipse. A sixth 
element is required to describe the satellite's position along the path at a 
particular point in time. [Ref. 2:p. 58] 

There are several sets of variables used to describe a satellites position. 
SMET uses the classical set of six orbital elements. The semi-major axis (a) and 
eccentricity (e) describe the size and shape of the ellipse. The inclination (i), 
longitude of the ascending node (h) and argument of perigee (g) describe the 
orientation of the ellipse are depicted in Figure 2-1. The epoch (to) is the time of 
perigee passage. It is used as the reference from which to compute the satellite's 
position with respect to time. 


B. OTHER STATIC PARAMETERS 


Based on the orbital elements, several other parameters are computed to help 
describe the orbit. These elements do not change with respect to time. The 
following equations show how to compute the semi-minor axis (b), period (T), 
perigee range (rp) and apogee range(ra). 


b 



T = 2ic'\/|' 

rp = a (1 - e) 
ra = a (1 + e) 


[Ref. 2:p. 33] 

[Ref. 2;p. 33] 

[Ref. 2:p. 25] 
[Ref. 2:p. 25] 
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Figure 2-1. Orbital Elements 
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A transformation matrix R is also computed using the longitude of the 
ascending node (h), argument of perigee (g) and inclination (i). This matrix is 
used in drawing the satellite and its orbital ellipse. It is also used to convert 
perifocal coordinates into geocentric-equatorial coordinates. 

cos(h)cos(g) - sin(h)sin(g)cos(i) -cos(h)sin(g) - sin(h)cos(g)cos(i) sin(h)sin(i) 
sin(h)cos(g) + cos(h)sin(g)cos(i) -sin(h)sin(g) + cos(h)cos(g)cos(i) -cos(h)sin(i) 

sin(g)sin(i) cos(g)sin(i) cos(i) 

[Ref. 2:p. 83] 


C. KEPLER PROBLEM 

Given a satellite's position with respect to its classical elements, it is fairly 
simple to compute the elapsed time since epoch. This is described at the end of 
Chapter III. It is, however, more difficult to compute a position given the 
elapsed time. First the mean anomaly (M) must be computed using the following 
equation. 


M 




(t - to) - 2k7t + Mo 


[Ref. 2:p. 220] 


SMET defines epoch (to) at perigee. Since the mean anomaly and true 
anomaly are both zero at perigee, the last term (Mo) will drop out from the 
previous equation. The next step, computing the eccentric anomaly (E), is more 
difficult. To do this, the following equation must be solved for E. It is, in 
classical terms, referred to as the Kepler problem [Ref. 2:p. 220]. 


M = E - e sin E 


[Ref. 2;p. 220] 
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This equation cannot be inverted in closed form to provide a function for E 
in terms of M. Therefore, another approach is required. The following Newton 
iteration scheme is used to solve for E. 

Eo = M 

Mn = En - e sin En [Ref. 2:p. 221] 

M - M 

r-e cos E„ [Ref.2:p.222) 

The second and third equations are repeated until M - Mn becomes 
sufficiently small. The true anomaly (0) can be then computed from E. 

® = [Ref.2:p.l87] 

The true anomaly is depicted in Figure 2-1. 

D. OTHER CHANGING PARAMETERS 

Once the position of the satellite (i.e., true or eccentric anomaly) is 
computed, other descriptions of its position can be computed. First the range (r) 
of the satellite is computed using the eccentric anomaly: 

r = a(l-ecosE) [Ref. 3:p. C-7] 

Next, the coordinates in the perifocal coordinate system are computed: 

P = r cos 0, [Ref. 2:p. 72] 

Q = r sin0, [Ref.2:p.72] 

W = 0 by definition. 
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The coordinates in the geocentric-equatorial coordinate system (X, Y and Z) 
are then computed with the transformation matrix R, 


x 


-p- 

Y 

-Z. 

= R 

Q 

LwJ 


[Ref. 2:p. 83] 


The last set of coordinates to be computed are the latitude and longitude or 
the sub-satellite point. Using simple geometry, we get 


lat = sin-1 



Ion = tan-i 

These coordinates are fixed in space. Note that the earth rotates at an 
angular velocity of co©. In order to compute the longitude with respect to the 
rotated earth, subtract the earth’s current rotation and allow for an arbitrary 
number of rotations to correct the previous equation. The corrected equation 
for longitude is 


Ion = tan-i 



CO© t + 2kK. 


Although the magnitude of the velocity (v) is not a description of the 
position, it does provide some insight as to the motion of the satellite. It is also 
instrumental in computing orbital maneuvers. The following equation provides 
the instantaneous velocity of the satellite. 



[Ref.3:p. C-11] 
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E. PERTURBATIONS 


A perturbation is a variation in the path of the satellite not explained by the 
Keplerian model for two-body motion used up to this point. Perturbations result 
from the gravitational influence of other bodies, atmospheric drag and the oblate 
shape of the earth. The two principle effects are the regression of the line of 
nodes and the rotation of the line of apsides [Ref 2:p. 156]. These two effects 
change the longitude of the ascending node and argument of perigee, 
respectively. When turned on, SMET uses the following equations to model 
these two principle perturbations. 


h = 


-3 r| 

2 a3-5 (1 - elf 


[Ref. 3:p. 4-4] 


• _ Vii Ji R© 
a3-5 (1 . e3)2 



[Ref. 3:p. 4-5] 


J 2 is the second order coefficient for perturbation accelerations due to the 
earth’s nonsphericity. Other perturbations are more computationally in ensive 
and provide changes to small to be visible with SMET. Therefore, SMET makes 
no attempt to model them. 


F. FIELD OF VIEW CALCULATION 

The Satellite View of SMET demonstrates how the altitude of a satellite 
effects the amount of the earth's surface visible to the satellite. This amount is 
referred to as the field of view (FOV) of the satellite. Note that the angle a in 
Figure 2-2 represents half of the FOV. If the range (r) is expressed in terms of 
earth radii (i.e., canonical units), the problem is simplified to: 

1 

sin a = - . 
r 
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Thus 
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III. SATELLITE LAUNCH 


A. GENERAL DESCRIPTION 

When a satellite is launched, it is normally attached to the top of a rocket. 
While the rocket is burning, the mass and velocity are constantly changing. 
Additionally, the effects of gravity and atmospheric drag change continuously. 
At the end of the rocket bum, the satellite will separate from the final stage of its 
rocket. The end effect of the rocket burn is to place the satellite at a particular 
point in space with a particular velocity. At this burnout point, the velocity 
vector of the satellite corresponds to a unique orbital path. SMET allows the 
user to enter parameters associatied this burnout point to start a satellite's 
motion. Modeling the complex interaction of the large number of variables that 
define the actual bum sequence is beyond the scope of SMET. 

B. PARAMETERS 

The values of the following variables at the time of burnout are used by 
SMET to compute the orbital parameters of the satellite. This is by no means the 
only way to define the burnout of the rocket. They were chosen because they are 
easy to visualize. 

Vo Velocity (with respect to a rotating earth) 

({) Latitude in earth-fixed coordinates 
X Longitude in earth-fixed coordinates 

Yo Flight Path Angle (measured upward from the local horizon plane) 
\|io Launch Azimuth (measured clockwise from North) 
alt Altitude above the surface of the earth 
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C. EARTH SPIN CORRECTION 

The earth rotates about its axis in 23 hours and 56 minutes. This produces a 
surface velocity at the equator of approximately 465.10 m/sec (1524 ft/sec) 
[Ref. 2;p. 306]. The surface velocity is a function of latitude, decreasing as the 
poles are approached. This explains why eastward launches or equatorial launch 
points are more fuel efficient. SMET adds this induced velocity to the satellite 
velocity defined by the burnout parameters. 

The velocity is broken up into components in a satellite fixed coordinate 
system. This is depicted in Figure 3-1. This system has a radial component, a 
North/South component and an East/West component. Using simple geometry, 
the following equations provide the solution for the velocity components. 

vr = Vo sin Yo 
vnvs = Vo cos Yo cos Xj/o 
Ve/Wo = Vo cos Yo sin Vo 
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The earth rotation component (vspin) is then computed from the latitude. R© 
is the radius of the earth and CO© is the angular velocity of the earth. 

Vspin = R©co© cos (}) [Ref. 2;p. 307] 

SMET adds vspin to ve/Wq and computes a new horizontal velocity 

component. The following equations, based on Figure 3-2, show the corrected 
values of velocity magnitude, flight path angle and azimuth. 


Ve/W = VEA^'o ''SPIN 





D. ORBIT PARAMETER CALCULATIONS 

SMET uses the corrected burnout parameters to compute the orbital 
parameters. The inclination is a function of the azimuth and latitude. 


i = cos'^ (sin \|/ cos (})) 

The range is a function of the altitude. 

The semi-major axis is a function of the range and velocity, 
r 


[Ref. 2:p. 142] 


a = 


0 . 


r 


[Ref. 3:p. C-5] 


The semi-latus rectum is a function of the range and horizontal velocity. 

P = [Ref. 3:p. C-71 

The eccentricity is a function of the semi-major axis and semi-latus rectum. 




[Ref. 3;p. C-6] 


The true anomaly is a function of the range and semi-latus rectum. 


e = u.n-.(^ 


[Ref. 3:p. C-9] 


The arguments of latitude and longitude are functions of the latitude and 
azimuth. 


sin 0 


u = tan-i . 

(cos <p cos \|f^ 

, /"sin (1) sin \ir 

(0 = tan-i -^^ 

1 cos y 


[Ref. 3:p. A-7] 
[Ref. 3:p. A-7] 
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The argument of perigee is a function of the argument of latitude and true 
anomaly. 


g = u-0 [Ref. 3:p. C-9] 

To compute the longitude of the ascending node, note that ^ - CO represents 
the distance between satellite's position in a rotational coordinate system and its 
position fixed coordinate system. Also note that (0©t represents the distance 
between an ascending node of 0° (fixed in space) and an ascending node that has 
rotated as much as the earth. 


h = - (0 + CO© t 

As discussed in Chapter II, SMET computes the true anomaly as a function 
of the elapsed time since epoch (to). Therefore, after the true anomaly has been 
computed, the eccentric anomaly (E) and mean anomaly (M) are computed using 
the eccentricity (e) and true anomaly. 


E = cos'^ 


^ e + cos 6 ^ 
^1 + e cos 0 J 


M = E - e sin E 


The last parameter, the epoch, is then computed using the mean anomoly, the 
semi-major axis (a) and current time (t). 
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IV. SATELLITE MANEUVERS 


A. GENERAL DESCRIPTION 

A satellite maneuver is executed by changing the velocity of the satellite. 
This change to the velocity vector can be in its magnitude, direction or both. A 
maneuver results in a new set of orbital parameters. A velocity change (Av) is 
modeled as an instantaneous impulse. Therefore, the satellite's position or range 
does not change while the velocity does. The Av vector can be broken up into 
component vectors (Av^, Avy and Ava) which cause different effects 
simultaneously. These three vectors cause a flight path angle change, a heading 
angle change and a change in the magnitude of the velocity vector, respectively. 
Both A\y and Avy change the direction of the velocity vector but do not change 
its magnitude. SMET computes the effects of each of the maneuvers 
independently and then updates the time of epoch. The time of epoch is 
computed such that the elapsed time will place the satellite in the same position 
based on the newly computed parameters. 

B. FLIGHT PATH ANGLE CHANGE 

A Avy will cause a change in the flight path angle (y). A flight path angle 
change is a maneuver in the orbital plane. The flight path angle is measured 
upward from the local horizon plane to the velocity vector [Ref. 3:p. 2-22]. The 
local horizon plane passes through the satellite and is perpendicular to the line 
passing through the center of the earth and the satellite. Avy is chosen such that 
the magnitude of the velocity vector is unchanged (i.e., vi = V 2 ). These vectors 
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form an isosceles triangle as can be seen in Figure 4-1. Given Ay, values for the 
eccentricity (e), argument of perigee (g), and true anomaly (0) can be computed. 



Figure 4-1. Flight Path Angle Change Vectors 


The eccentricity changes as the impulse changes the shape of the orbit. The 
argument of perigee and true anomaly change as a result of the rotation of the 
line of apsides. This can be seen in Figure 4-2. Since A\y has no out-of-plane 
component, there is also no change in the inclination (i) or the longitude of the 
ascending node (h). The following equation shows that the semi-major axis (a) 
does not change since the range and the magnitude of the velocity do not change. 

a =-^ [Ref.3:p.C-5] 

2 - 

The eccentricity can be computed given Ay by combining the following two 
equations for the cosine of the flight path: 


vnp 

cosy =—“, 
' rv 


[Ref. 3:p. 3-15] 


cos yi = cos yi (cos Ay - tan yi sin Ay). 
This results in 


VupT _ 

r2V2 rivi 


(cos Ay - tan yi sin Ay). 


[Ref. 3:p. 3-15] 















Recalling that the range and magnitude of the velocity are unchanged, they 
can be removed from the previous equation. 

= (cos Ay - tan 7 i sin Ay) 

The semi-latus rectum is computed with 

p = a (1 - e2). [Ref. 2:p. 24] 

Substituting this into the previous equation while recalling that the semi¬ 
major axis remains unchanged, results in 

= (cos Ay - tan 71 sin A 7 ), 

1 - ef = (1 - ei)(cos Ay - tan 71 sin Ay)^. 

The tangent of the flight path angle can be computed from the eccentricity 
and eccentric anomaly. 

e, sin E, 

tan 7 i =- ] ...... ’ [Ref. 3:p. C-10] 

Vl - e? 

Thus a new eccentricity can be computed with Ay, ej and Ej. 



e 2 = 


V 


1 - (1 - f^) 


cos Ay 


0 ] sin E| sin Ayy^ 




e 2 

Cl 




Another equation relating the flight path angle to the satellite's position 
provides the means for computing the new true anomaly. 


0 = tan 


i^pjan^ 


[Ref. 3:p. C-9] 

Using tan 7 i, which was computed earlier, we can compute the tangent of the 
new flight path angle. 


tan 7 , -I- tan Ay 
tan 7 , = ^ 

1 - tan 7 i tan Ay 


[Ref. 3:p. 3-19] 
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Using P 2 = a (1 - 62 ), the new anomaly can thus be computed. 



Since the satellite has not moved and the line of nodes has not changed, the 
argument of latitude (u) remains unchanged. Knowing that g = u - 0, we can 
compute the new argument of perigee using the old argument of perigee and old 
true anomaly with the new true anomaly [Ref 3:p. 3-8], 


g2 = gi + 01 - 02 

C. HEADING ANGLE CHANGE 

A Avy will cause a change in the heading angle (y). A heading angle change 
is a maneuver out of the orbital plane. The heading angle is the angle in the local 
horizon plane measured from the line through the satellite and the north point 
and the velocity vector [Ref. 3:p. 2-22]. Again, Avy is chosen such that the 
magnitude of the velocity vector is unchanged (i.e., vi = V 2 ). These vectors 
form an isosceles triangle as can be seen in Figure 4-3. With Ay, values for the 
inclination (i), argument of perigee (g) and longitude of the ascending node (h) 
can be computed. 
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The inclination changes as the orbital plane rotates about the line connecting 
the center of the earth and the current satellite position. This rotation can be 
seen in Figure 4-4. This rotation also changes the point where the orbital plane 
intersects with the equatorial plane, causing the ascending node to change. 
Finally, as the plane rotates, the distance between the ascending node and the 
current position as measured along the orbital path (the argument of latitude (u)) 
changes, causing the argument of perigee to change. Since A\}r has no in-plane 
component, there is no change to the shape or in-plane orientation of the orbit. 
Thus, the eccentricity (e), semi-major axis (a) and true anomaly (0) remain 
unchanged. 

Given A\|/, the new inclination and argument of latitude can be computed 
using the following equations: 


cos h = cos ii cos Ay -t- sin ii cos Uj sin Ay, [Ref. 3:p. 3-6] 
sin 12 cos U 2 = sin ii cos Ui cos Ay - cos ii sin Ay, [Ref. 3:p. 3-6] 
ui = gi + 01. [Ref. 3:p. 3-8] 

The new inclination and argument of latitude are found by rearranging the 
previous equations. 


i 2 = cos-1 (cos ii cos Ay + sin ii cos(gi+0i) sin Ay ) 
(sin ii cos(gi+0i) cos Ay - cos ii sin Ay 

sin i 2 


e 


■) 


U2 = COS-1 

Since the true anomaly is unchanged (i.e., 0i = © 2 ), the new argument of 
perigee is simply; 


g2 = U2-0i. 
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Figure 4-4. Heading Angle Change 
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The new longitude of the ascending node is computed by noting that the 
change in the longitude of the ascending node is equal to the change in the 
argument of longitude [Ref. 3:p. 3-11]. The argument of longitude is computed 
as follows: 


CO = tan-^ 


/sin u cos h 
cos u / 


[Ref. 3:p. A-8] 


™ , /^sin ui cos iA , ^in U2 cos i2^ 

Thush2 =hi -i-tan-i ——7--^j-tan-i —-±- - 

^ cos Ui y cos U2 / 


D. VELOCITY MAGNITUDE CHANGE 

A Avq will change the magnitude of the velocity (v 2 = vi + Avq) without 
changing the direction (Figure 4-5). Based on the new velocity magnitude, 
values for the semi-major axis (a) and eccentricity (e) can be computed. These 
changes are reflected as a change in the shape of the orbit (Figure 4-6). Since 
the heading angle and flight path angle are unchanged, the parameters defining 
the orientation of the orbital plane (i.e., i and h) are unchanged [Ref 2.:p. 3-20]. 
If the orbit was circular prior to the change, the argument of perigee (g) and 
true anomaly (6) must also be set, since those two parameters have little meaning 
in a circular orbit. The velocity change at a position other than perigee or 
apogee will also cause a rotation of the line of apsides. Therefore, the argument 
of perigee and true anomaly must be computed as they were for a flight path 
angle change. 



Figure 4-5. Velocity Magnitude Change Vectors 


23 








Figure 4-6. Velocity Magnitude Change 
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The semi-major axis and the eccentricity are computed with the following 
equations: 

[Ref. 3:p. 3-20] 


32 = 


r v^’ 




62 




2 

Vh2 


[I 32 


where Vn 2 = V 2 cos y. [Ref. 3:p. 3-20] 


The flight path angle (y) does not change since the velocity vector was not 
rotated. Therefore, the old values of e and E can be used to compute Y- 


e, sin El 
tany =—'-^ 


[Ref. 3:p. C-10] 


Vl - ei 

If the original orbit was circular (i.e., ei = 0), the fact that the argument of 
latitude remains unchanged is used to compute a value for the argument of 
perigee and true anomaly. The current position will be perigee if there was an 
increase in velocity. 


g 2 = u = gi - 1-01 
02 = u - g2 = 0° 

If the velocity decreased, the satellite will be at apogee of the new orbit. 

g 2 = u- 180°= gi -t- 01 - 180° 

02 = U - g2 = 180° 

E. UPDATING EPOCH 

Changes in the magnitude of the velocity and flight path angle cause a change 
to the true anomaly (0). As in the satellite launch equations, the jpoch is 
computed as a function of the true anomaly and the current time with respect to 
the eccentricity and semi-major axis of the new orbit. 
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F. IMPULSE VECTOR DIVISION 

Given a velocity change vector (Av) with an associated yaw and pitch, the in¬ 
plane and out-of-plane components (Avp and Avop) can be computed. Using Avp 
and Avop, the desired vector components (Avy, Avy and Ava) can be computed. 
Unfortunately, Avy, Avy and Aw a are not orthogonal vector and, therefore, are 
not easily computed. Figure 4-7 depicts the relationship between the two sets of 
vectors. The same set of vectors is depicted again as part of Figure 4-10. 

Earlier sections noted that A\y and Av,,, formed isosceles triangles when 
added to the velocity vector. These isosceles triangles had associated angles (Ay 
and A\|/) which were used to compute the flight path change maneuver or heading 
change maneuver. Similarly the new velocity (V 2 ), and not Ava, is actually used 
to compute the new orbital parameters for a velocity magnitude change 
maneuver. SMET, therefore, does not use Avy, Av^ and Avq, SMET actually 
computes Ay, Ayandva. 




The yaw associated with a maneuver is used to compute the in-plane and out- 
of-plane components of the Av vector. 

Avp = Av cos(Y aw) 

Avop = Av sin(Yaw) 

The pitch associated with the maneuver is used to compute that portion of the 
in-plane component (Avp) which is in the direction of motion and that portion 
which is perpendicular to the direction of motion. These vectors are added to 
the original velocity as depicted in Figure 4-8. 



Figure 4-<S. Computing Flight Path Angle Change 


This right triangle provides the means with which to compute the change to 
flight path angle (Ay). 

vi + Avp cos(Pitch) J 

The hypotenuse of the right triangle is computed for later use. 

v'l + Ava = V(Avp sin(Pitch))2 -h (vi + Avpcos(Pitch))2 
In order to compute the heading angle change (Ay), the out-of-plane 
components of Av (Avop) are added to the previously computed in-plane 
components (v'l + Ava). Figure 4-9 shows that a right triangle is formed again. 
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Figure 4-9. Computing Heading Angle Change 


This right triangle provides the means with which to compute the change to 
the heading angle 

AY = tan-i f - V - ? . 1 
l^vi + Avaj 

The hypotenuse of this right triangle (v’i + Ava) is V 2 . This can be seen in 
Figure 4-10 which combines the two previous figures. 

V2 = vi + Ava = VAvop^ + (v’l + AvaP 
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V. USER’S GUIDE 


A. OVERVIEW 

This chapter describe SMETs display and user interface. Figure 5-1 shows 
the five windows used by SMET - the Main Window, two Auxiliary Windows, 
the Parameter Window and the Satellite List Window. The Main Window is a 
large square window at the center of the screen. There are two smaller 
Auxiliaiy' Windows below the Main Window. These three windows are called 
viewing windows. The Parameter Window is a long narrow window on the left 
side of the screen. The Satellite List Window is a long narrow window on the 
right side of the screen. 

SMET updates all satellite positions continuously with respect to time. One 
satellite is designated as the Current Satellite. This satellite is listed at the top of 
the Satellite List Window and has its parameters displayed in the Parameter 
Window. Additionally the satellite viewpoint, discussed later, belongs to the 
Current Satellite. All user inputs affect the Current Satellite. Any satellite can 
become the Current Satellite with the use of the Satellite List Window. 

B. VIEWING WINDOWS 

The three viewpoints discussed in this section can appear in any of the three 
viewing windows. These images are completely interchangeable. The user can 
switch an image from an Auxiliary Window into the Main Window by simply 
pressing the middle mouse button while the cursor is in the desired Auxiliary 
Window. The image that was in the Main Window will be moved to the selected 
Auxiliary Window. 
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1. OmniViewer View 

The OmniViewer View allows the user to observe all satellites as they 
orbit around a rotating earth model as seen in Figure 5-2. Axes for the 
geocentric-equatorial coordinate system are displayed for user orientation. An 
orbital path is displayed for each satellite. The user can choose to draw it as a 
colored line or as a colored transparent ellipse. Normal and perigee vectors can 
also be displayed for any satellite to assist in user orientation. 

The user's viewing position is referred to as the Viewer X, Y and Z in 
the Parameter Window. Changing the user’s viewing position is discussed in 
SMET System Parameters. It can be set automatically or adjusted by the user. 
When Automatic OmniViewer is turned on, the viewer will be positioned far 
enough away to be able to see all orbital paths. The viewer position will be 
adjusted automatically whenever a satellite is added or any orbital path is 
changed as a result of a maneuver or user input. During satellite maneuvers or 
launches, Automatic OmniViewer can be turned off. This allows the user to see 
the event on a single scale. 

2. Satellite View 

The Satellite View allows the user to see that portion of the earth visible 
to the Current Satellite. The image is a "fish-eye lense" type of view. The user 
sees a circular image of the earth encircled by a ring colored in the Current 
Satellite's color. The center of the image corresponds with the Current 
Satellite's sub-satellite point. The outer edge represents the horizon as seen from 
the altitude of the Current Satellite. As the Current Satellite’s altitude changes, 
the field of view is adjusted to maximize the size of the image. Therefore, the 
size of the image does not generally change. Figure 5-3 shows a typical view 
from a satellite in a Molniya orbit. 
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Figure 5-2. OmniViewer View 
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Figure 5-3. Satellite View 
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3. Mercator View 

The Mercator View shows a map of the world as seen in Figure 5-4. 
The sun’s relative plot is depicted as a yellow transparent ellipse on the map. 
This ellipse does not attempt to represent day/night transitions. Each satellite 
will have its sub-satellite point plotted in its color every time the satellite’s 
position is updated. Satellite Parameters discusses how to turn plotting on. 

C. PARAMETER WINDOW 

The Parameter Window displays names and values of parameters for the 
Current Satellite and the SMET system. The user can change the value of many 
parameters with mouse and keyboard inputs. Distance and velocity parameters 
are displayed in one of three types of units. The three types are canonical (DU, 
DU/TU), metric (km, km/sec) and english (NM, ft/sec). 

1. Satellite Parameters 

There are five sets of Satellite Parameters - Classical, Secondary, 
Position, Options and Combined. The Classical set includes the parameters that 
describe the shape and orientation of the orbit. Most of these parameters can be 
changed by the user. The Secondary set also includes parameters that describe 
the shape of the orbit. Most of these parameters are actually computed by SMET 
for display and cannot be changed by the user. The Position set includes 
parameters that describe the satellites position in various coordinate systems. 
These parameters are also computed and cannot be changed by the user. The 
Options set is a set of parameters that SMET uses. 'The Combined set includes all 
the members of the ether sets listed above. All parameter sets display the 
Current Satellite’s name at the top. Figure 5-5 shows the Combined set of 
satellite parameters with the system parameters listed at the bottom. 
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Figure 5-4. Mercator View 
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Satellite Parameters: 


Name MoUiiva Orbit 
Semi-Major 56561.78 Ic 
Semi-Minor 17568.97 k 
Eocentrictty075 
Inclination 63.43 
Ascending Node 99.55 
Arg of Periapsis 270 00 
True AnomaK’ 180.97 
Epoch 25;54:(k) 

Period 11:58.02 
Apogee 4648312 km 
Perigee 6640,45 km 
Range 46463.18 km 
Velocitv’ 147 km/sec 
X 18308 79 km 
Y 9851.25 km 
Z 41551.97 km 
P -46456.53 km 
Q -786.01 km 
Latitude 63.42 N 
I,ongjt'i<fe 15009 W 
Color Red 
Fill Off 
Vectors On 
Plot On 


System Parameters: 
Clock On 

Current Time 44:30:00 
Dock Speed 00:30:00 




Political B«diefs On 
PciturlMrtioiis Off 
Auto Omni Viewer On 
Viewer X 46483.12 km 
Y 637815 km 
Z 6378.15 km 
Units Metric 


Figure 5-5. Parameter Window 
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a. Classical Set 

The Classical set of satellite parameters includes the semi-major 
axis, semi-minor axis, eccentricity, inclination, longitude of the ascending node, 
argument of periapsis, true anomaly, and time of epoch. The true anomaly is 
computed by SMET as a function of the elapsed time since epoch. It cannot be 
changed by the user. All other parameters can be adjusted by the user as 
described below. 

b. Secondary Set 

The Secondary set of satellite parameters includes the perigee 
distance, apogee distance, time of epoch, orbital period, current range and 
velocity. Except for the time of epoch, all of these parameters are computed by 
SMET and cannot be changed by the user. 

c. Position Set 

The Position set of satellite parameters displays the current position 
of the satellite in three different sets of coordinate systems. The X, Y and Z 
coordinates correspond with the geocentric-equatorial coordinate system 
displayed in the OmniViewer window. The P and Q coordinates correspond 
with the perifocal (P-Q-W) coordinate system. Both of these sets of coordinates 
are displayed in the current unit type (ex. metric). The Latitude and Longitude 
coordinates correspond with the earth-fixed coordinate system seen in the 
Mercator View window. None of these parameters can be change by the user. 

d. Options Set 

The Options set include the color of the satellite and its ellipse, if 
the orbital ellipse should be filled or empty, if the normal and perigee vectors 
should be displayed and if the satellite should be plotted on the Mercator map. 
All of these parameters can be changed by the user. 
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e. Combined Set 

The Combined set includes all the members of the other sets listed 

above. 

2. SMET System Parameters 

SMET system parameters are listed below the set of satellite 
parameters. This set of parameters includes the clock, current time, clock speed, 
detailed maps, political borders, perturbations. Automatic OmniViewer, 
OmniViewer X, Y, and Z, and units. The clock parameter is used to start and 
stop the clock. Current time displays the elapsed time in hours:minutes:seconds 
format. It uses 00:00:00 GMT on the day of Vernal Equinox as the reference 
point. Clock speed allows the user to set the number of minutes between position 
updates. It is initially set to 30 minutes. The detailed maps parameter allows the 
user to switch between detailed and coarse globe and map data. The political 
borders parameter allows the user to turn off the display of political borders on 
the globe and map. Both of these options allow a user with a large number of 
satellites to speed up SMET by drawing less detailed images. They are both 
initially set to on. The perturbations parameter turns on the regression of the 
line of nodes and the rotation of the line of apsides. This is initially set to off. 
Automatic OmniViewer positioning can be turned on and off. It is initially set to 
on. If any of the OmniViewer coordinates (X, Y or Z) are changed. Automatic 
OmniViewer positioning is turned off. The OmniViewer coordinates are 
geocentric-equatorial coordinates. They are displayed in the current unit type. 
The units parameter displays the current unit type (canonical, metric or english). 
It is initially set to canonical. All system parameters can be changed by the user. 
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3. Parameter Input 

A parameter in the Parameter Window is selected by placing the cursor 
over it and pressing the middle mouse button. If it can be changed, a dialog box 
or menu will appear. Dialog boxes are answered with keyboard entries. All 
keyboard entries must be entered in the current unit type (e.g., metric). A menu 
is answered by selecting an entry and pressing the right button. Parameters that 
toggle between "On" and "Off" will change with no other input required. 
Similarly, the units parameter will cycle through "Canonical", "Metric" and 
"English" each time it is selected. Certain parameters are computed by SMET. 
If one of these parameters is selected, the user is told that it cannot be changed. 

SMET also allows the user to adjust certain parameters with five dials 
on the dial box. These include the clock speed and four of the classical 
parameters. Figure 5-6 shows how the dials are arranged. The dials are used 
for user convenience and simple demonstrations of the manipulation of 
parameters. If no dial box is available, these variables can still be input by using 
the Parameter Window. Whenever a dial is turned, the appropriate parameter 
will be changed. The clock speed dial is initially calibrated to operate between 
zero and six hours. If a clock speed greater than six hours is desired, the user 
can enter it by using the Parameter Window. 

A change to a satellite's classical elements will radically alter its 
position. Therefore, SMET erases all displayed plots for the Current Satellite 
and begins plotting again. Similarly, if the current time is changed, all plots for 
all satellites are erased. 
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D. SATELLITE LIST WINDOW 

The Satellite List Window, seen in Figure 5-7, lists all of the satellites in the 
system. The Current Satellite is listed at the top of the list. The user changes the 
Current Satellite by placing the cursor over the desired satellite and pressing the 
middle mouse button. The other satellites will shift down one position as 
necessary to reorganize the list. 

E. MENU COMMANDS 

When the user presses the right hand button on the mouse, a menu will 
appear near the position of the cursor. A menu item is selected by pressing the 
right hand mouse button again while the cursor is over a particular command or 
option. Certain menu items will have an arrow to the right. These items 
indicate that another menu (a rollover menu) will appear if the cursor is moved 
to the right of that menu item. Pressing the right hand mouse button while it is 
over the menu title or at any time that it is not over a menu item will abort the 
process of making a menu selection. 

The main menu allows the user to add a satellite, maneuver the current 
satellite, delete the Current Satellite, save satellites' parameters into a file, take a 
picture, choose the set of satellite parameters to display in the Parameter 
Window, list the user instructions or exit the program. 

1. Add Satellite 

The Add Satellite command has three options that appear in a rollover 
menu. The user can define an orbit's parameters, define the parameters for a 
launch, or read a file of previously saved satellites. 
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SMET 


Major Eccentricity 


Inclination Ascending 

Node 


Argument Clock 

of Perigee Speed 


Figure 5-6. Dial Box 
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SMET Satellite List: 

Molniya OilMt 
LowLarth O^t 
Ge<»tationaiy Orbit 
Gwsynchronoys Orbit 
Another Orbit 
Yet Another Orbit 


Figure 5-7. Satellite List Window 
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a. Define Orbit 

Defining an orbit's parameters also has a rollover menu. Here the 
user is given the option of selecting from three pre-defined orbits (a Low Earth 
orbit (300 NM), a Molniya orbit, and a Geosynchronous orbit). The user also 
can select "Input Parameters". This allows the user to enter the classical 
parameters of the desired orbit through a series of dialog boxes. Parameters 
must be entered in the current unit type (e.g. km, if in metric). The user is 
asked to enter a name for the new satellite in all cases. It will be given a 
randomly selected color and have plotting, fill, and vectors turned off. The new 
satellite becomes the Current Satellite. 

b. Define Launch 

The Define Launch command also initiates a series of dialog boxes. 
The user is asked to enter a name for the new satellite and burnout parameters 
related to the launch. These parameters include the time at which burnout occurs 
(not tc be confused with the length of the rocket bum), the altitude above the 
surface of the earth (in the current unit type, e.g. km), the velocity (in the 
current unit type, e.g. km/sec), the earth-fixed burnout latitude and longitude 
(±degrees), the flight path angle and the launch azimuth. If the time of burnout 
is earlier than the current time, the satellite will be placed at a position where it 
would have been had the launch been executed at the appropriate time. If not, a 
launch record is saved until the appropriate time of execution. When executing a 
launch, SMET will add in the appropriate earth rotation correction. The new 
satellite will be given a randomly selected color and have plotting, fill, and 
vectors turned off. It will become the Current Satellite. See Chapter HI for 
more information. 
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c. Read File 

The Read File command will read a satellite file previously saved 
using the Save Satellites command discussed below, 

2. Maneuver Satellite 

The Maneuver Satellite command initiates a sequence of dialog boxes. 
The user is asked to enter parameters for a maneuver of the Current Satellite, 
The first parameter is the true anomaly at which the maneuver should occur. If 
there are any maneuvers pending for that satellite, the user will be asked to enter 
the number of degrees of travel following a previous maneuver. When a 
previous maneuver is executed, this quantity will be added to the true anomaly at 
the end of the previous maneuver. For example, if 900 is entered the new 
maneuver will occur two and a half revolutions after the previous maneuver 
occurs. There is no limit to the number of maneuvers which can be added to a 
satellite's queue. Maneuvers will be executed in the order in which they are 
entered. A satellite does not have to be the Current Satellite at the time of the 
maneuver. 

After the first entry the user is asked to choose between two methods of 
entering the parameters of a maneuver. The first option is to enter the magnitude 
of the velocity change along with an associated pitch and yaw. The user enters 
the magnitude in the current units (e.g., km/sec). The pitch and yaw are entered 
in degrees. Pitch is a right-hand rotation about the satellite’s z-axis (thumb 
pointing towards the earth). Yaw is a right-hand rotation about the satellite’s 
y-axis (thumb pointing opposite to the orbital plane's normal vector). SMET 
uses the magnitude, yaw and pitch of the velocity change to compute the changes 
to the magnitude of the velocity, flight path angle and heading angle of the 
satellite at the time of the maneuver. The second option for entering a maneuver 
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is to enter the changes to the magnitude of the velocity, fl'ght path angle and 
heading angle directly. This option is useful to users who are more interested in 
the affects of a maneuver than in the total Av required. See Chapter IV for more 
information. 

3. Delete Satellite 

The Delete Satellite command will remove the Current Satellite from 
the system permanently. The second satellite listed in the Satellite List Window 
becomes the Current Satellite. 

4. Save Satellites 

The Save Satellites command will first ask the user for a file name. The 
user may enter any Unix path and file name. If the file name is not unique, the 
old file will be overwritten. The user will then be asked for a number of 
satellites to be saved. Satellites are saved from in the order that they appear in 
the Satellite List Window. If the user only desires a specific subset of the 
satellites in the system, each satellite desired should first be selected in sequence 
until they are at the top of the list. Any file saved using this procedure can be 
read by Read File command discussed above. Saving satellites does not remove 
them from the system. 

5. Take Picture 

The Take Picture command has four options that will appear in a 
rollover menu. All options request a file name. The first two options will save 
an sgi-rgb formatted file. These files can be shown and manipulated with several 
image tools available on the IRIS. These image tools include showsgi, shrinksgi 
and others. They are available in the imagctool subdirectory of SMETs 
directory. The first option saves the entire screen as it appears. The second 
saves the entire screen after the background has been changed to white and white 
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items have been changed to gray. The third and fourth options will save the 
entire screen or the Main Window in Encapsulated PostScript (EPS) format. 
These files can be printed on a PostScript laserwriter or read by any software 
capable of reading EPS images. SMET will remain frozen for 30 to 90 seconds 
while the image is being saved. 

6. Satellite Parameters 

The Satellite Parameters command has a rollover menu. This menu 
allows the user to select which of the five sets of satellite parameters will be 
displayc ' in the Parameter Window. 

7. Instructions 

The Instructions command will temporarily halt SMET and display the 
user instructions. 

8. Exit 

The Exit command will quit SMET. The user is offered an opportunity 
to save any sateLites desired. Any satellites in the system that have not been 
saved w'ill be lost. 
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VI. PROGRAM MODIFICATIONS 


A. GENERAL 

There are several dummy function calls entered into SMETs source code. 
These dummy procedures can be used to make minor modifications to SMET 
without making major changes to the code. This capability is directed primarily 
at thesis students requiring a specific simulation. All of the functions described 
below are included in the file smet_user.c which is included in Appendix C. This 
file can be modified by novice C programmers for desired modifications to 
SMET. Other smet_xxx.c files should not be changed by inexperienced 
programmers. Modified versions of SMET should have the name altered to 
prevent confusion with the original SMET. 

B. USER PROGRAM INITIALIZATION ROUTINE 

The procedure userjnit can be used to read in satellite files and initialize 
program variables to values different than the default value. 

C. USER SATELLITE INITIALIZATION ROUTINE 

The procedure user_new_sat can be used to modify any satellite variables 
to values different than the default value. The SMET satellite structure provides 
four unused variables (user_floatl, user_float2, user_intl and user_int2) which 
can be used in a modified program for such things as orbital perturbations or 
ASAT simulations. 
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D. USER PROGRAM UPDATE ROUTINES 

SMET provides two procedures that can be used to check for an exit 
condition or a condition for executing a launch or maneuver. Procedures 
user_pre_update_satellites and user_post_update_satellites can also be 
used to check the elapsed time or to compare the positions of the satellites. One 
routine is called before all of the satellites are updated and the other is called 
after. 

E. USER SATELLITE UPDATE ROUTINES 

SMET provides two procedures that can be used to perturb an orbit or 
delete a satellite based on a given set of parameters. The procedure 
user_pre_orbital_parameters is called before the satellite's changing orbital 
parameters are updated. The procedure user_post__orbital_parameters is 
called after the changing parameters have been updated. Both procedures take a 
satellite and the time as arguments. 
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VII. SUMMARY, RECOMMENDATIONS AND FUTURE 

DIRECTIONS 


SMET provides beginning students with the ability to learn the difficult 
concepts of orbital mechanics with greater ease. Satellite maneuvers, which are 
not covered in introductory courses in great detail, can now be taught with 
simple presentations. Additionally, a means of generating simulations is 
available for more advanced students. 

In order to achieve the maximum benefit of this project, the Space Systems 
Academic Group should purchase an IRIS graphics workstation capable of 
running SMET. Several students in the Space Systems Operations Curriculum 
have already received instruction in programming the IRIS. Already, one 
student is creating a smet_user module for a SMET based ASAT simulation. As 
a start, modules could be written for standard transfer orbit maneuvers or 
perturbations not currently modeled. Space Systems Engineering students could 
also benefit from the use of a CAD system or other engineering software 
available for the IRIS. As an interim measure, the Space Systems Academic 
Group should negotiate limited access to an IRIS belonging to another 
department for the development of class presentations. 

It is hoped that the software will be used elsewhere and continue to evolve. 
Future enhancements not currently provided include the ability to read 
NAVSPASUR Charlie Element formatted files for entering satellites, an 
illumination model to represent the sun, a more accurate representation of 
perturbations and a capability to automatically compute orbit transfer or 
intercept requirements. 
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APPENDIX A - COMMON VARIABLES 


Variable 

Name 

Formula 

Name 

Program 

Name 


Semi-Major 

Axis 

a 

major 

Description of conic section 


e 

ecc 

Description of conic section 

Semi-Latus 

Rectum 

p 

P 

Description of conic section 

Range 

r 

range 

Distance from the center of the earth to the 
satellite 

Velocity 

V 

velocity 

Magnitude of the velocity of the satellite 

Inclination 

i 

incl 

Angle between the plane of the equator and 
the plane of the orbit 

Longitude of 
the Ascending 
Node 

h 

asc 

Angle in the equatorial plane between the 
right ascension of the sun at Vernal Equinox 
and the right ascension of the satellite 


g 

peri 

Angle in the orbit plane from the line of 
nodes to perigee 

IPSHi 

0 

anom 

Angle in the orbit plane between perigee and 
the satellite 

Eccentric 

Anomaly 

E 

E 

Angle in the orbit plane between perigee and 
a "projection" of the satellite's position onto a 
circumscribed auxiliary circle 

HMI 

M 

M 

The product of the mean motion and elapsed 
tin" 

Epoch 

to 


Tiriic of perigee passage 

Heading 

Angle 

V 


Angle in the local horizon plane between the 
line through the satellite and the North point 
and the line of motion measured 

Flight Path 

Y 

f_path 

Angle in the orbit plane measured upward 
from the local horizon plane to the line of 
motion 

ri|n‘|iiinniit^ 

(0 


Angle in the equatorial plane between the line 
of nodes and the meridian of the satellite 

Argument of 
Latitude 

u 

u 

Angle in the orbit plane between the line of 
nodes and the satellite 

Latitude 


lat 

Latitude of the subsatellite point 

Longitude 

■■ 

Ion 

Earth-fixed Longitude of the subsatellite 
point 

Gravitational 

Parameter 



Constant for computation of earth orbits 
(1 DUffi3/TU®2 = 1.407647 x 1016 ft3/sec2) 
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APPENDIX B - VIDEO RECORDING AND PRINTING 


The following instructions apply to the implementation of SMET in the 
Graphics and Video Laboratory located in the Computer Science Department of 
the Naval Postgraduate School. 

Video Recording - With Professor Zyda's permission, an S-VHS or a 
standard VHS tape can be recorded while running SMET on IRISl. This tape 
can then be edited and dubbed as necessary for presentations. Note that some 
resolution will be lost because standard video signals have a lower resolution 
than the IRIS. The top and bottom edges of the image will also be clipped 
because of the different aspect ratios. The YEM scan converter should be turned 
on in order to record. The video recorder will record everything appearing on 
the video monitor. 

Printing EPS Files - EPS files are saved using one of the EPS options of 
the Take Picture command. These files can be printed on any PostScript laser 
printer. The file can be printed with the following command: "Ipr -Pps2 
filename". The file can also be imported into other programs such as 
Framemaker for incorporation into other documents as was done in the 
preparation of this document. 
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APPENDIX C - SMET PROGRAM LISTING 


SMET consists of the following files: 

smet.h (p. 54 - 59) contains all constants and structures used by 
SMET. Also has all forward function, procedure and global 
variable declarations. 

smet_vars.c (p. 60 - 64) defines all global variables. 

smet_addelete.c (p. 65 - 68) includes procedures to add, delete 
and switch satellites. Additionally, procedures to add and 
delete plots to satellite structures are included. 

sniet_axes.c (p. 69) includes procedure to draw the axes in the 
OmniView window. 

smet_dialog.c (p. 70 - 72) includes procedures for dialog and 
message boxes. 

smet_files.c (p. 73 - 77) includes procedures that read and write 
satellite files, as well as reading SMET data files and writing 
image files. 

smet_globe.c (p. 78 -79) includes procedures to draw the globe in 
the OmniView and Satellite View windows. 

smet_init.c (p. 80 - 82) includes SMET initialization procedures. 

smet_kepler.c (p. 83 - 85) includes procedures that update a 
satellite's orbital parameters. 

smet_launch.c (p. 86 - 87) includes procedures to launch a 
satellite. 


52 








sniet_main.c (p. 88 - 89) is the heart of SMET. It includes the 
main execution loop. 

smet_maneuver.c (p. 90 - 95) includes procedures to maneuver a 
satellite. 

smet_map.c (p. 96 -97) includes procedures to draw the map in 
the Mercator View window. 

smet_orbit.c (p. 98) includes procedure to draw a satellite's 
orbital path in the OmniView window. 

smet_param.c (p. 99 - 104) includes procedures that handle 
Parameter window events. 

smet_queue.c (p. 105 - 109) includes procedures read the queue 
and handle events not covered in other files. 

smct_sat.c (p. 110 - 113) includes procedures to draw a satellite 
in the OmniView window, 

smet_user.c (p. 114) includes dummy procedures discussed in 
Chapter VI. 

smet_windows.c (p. 115 - 123) includes procedures that draw the 
various windows. 

makefile (p. 124) can be used to compile SMET in the NPS 
Computer Science Department Graphics and Video 
Laboratory. Files in the imagetools sub-directory will also be 
required. 

smet_data.d (not included) is a binary data file necessary to draw 
the globe and map. 
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it.h Mon S«p 3 17:59:00 1990 1 


/**********************************/ 
{* */ 

/ * smet. h * / 

/» by Carlos I. Noriega */ 

/* */ 

/* This file contains the global variables, constants and */ 

/* structures used by the SMET demonstrator. */ 

/* */ 

/* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/* */ 

/****•******★*********************■**/ 

♦include <math.h> /* math functions file */ 

♦include <string.h> /* string routines */ 

♦include "gl.h" /* graphics library file */ 

♦include "device.h" /* device header files */ 

♦include "stdio.h" /* I/O routines */ 

♦include "fmclient.h" /* IRIS font manager */ 

♦include "image3upport/image_types.h" /* for picture saving */ 

♦include "imagesupport/image_funcs.h" /* for picture saving */ 


♦define FINEPOINTS 12843 
♦define COARSEPOINTS 3322 
♦define BORDERPOINTS 1407 

♦define OMNISATMERC 1 /* Which view is in which window? */ 
♦define OMNIMERCSAT 2 /* These constants used by switch_window3() */ 
♦define SATMERCOMNI 3 ~ 

♦define SATOMNIMERC 4 
♦define MERCOMNISAT 5 
♦define MERCSATOMNI 6 

♦define EARTH 8 /* Color Codes */ 

♦define COAST 9 
♦define BORDER 10 
♦define LATLONG 11 
♦define SUN 12 
♦define SAT_RED 13 
♦define SAT_GREEN 14 
♦define SAT~YELLOW 15 
♦define SAT_BLUE 16 
♦define SAT_MAGENTA 17 
♦define SAT_CYAN 18 
♦define SAT_WHITE 19 
♦define SAT_GPAY 20 
♦define MAX COLORS 21 


♦define CANONICAL 0 /* Unit Codes */ 

♦define METRIC 1 

♦define ENGLISH 2 


♦define DISTANCE 1 /* 


♦define SPEED 2 
♦define DEGREE 3 
♦define LAT 4 
♦define LON 5 
♦define TIME 6 
♦define FLAG 7 
♦define FLOAT 8 
♦define STRING 9 


♦define COLORCODE 10 
♦define UNITCODE 11 


for parameter type field */ 
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SMt.h 

Mon Sep 3 

17:59:00 1990 

2 




♦define 

ROTATION_RATE 

4.37526951286e-3 

/* 

earth rate 

in rads/min 


♦define 

TU 

13.44686457 

/* 

minutes per TU 

*/ 

♦define 

DU_METRIC 

6378.1459997 

/* 

KM per DU 


*/ 

♦define 

DU ENGLISH 

3443.9227864 

/* 

NM per DU 


*/ 

♦define 

DU_PER_TU_METR1C 

7.9053682830 

/* 

KM/SEC per 

DU/TU 

*/ 

♦define 

DU_P ER_TU_ENGL1SH 

25936.24764 

/* 

FT/SEC per 

DU/TU 

*/ 


♦define 

PI 

3.14159265359 

/* 

Not mom's apple ... 

*/ 

♦define 

TO_DEG 

57.2957795132 

/* 

Convert degrees to radians 

*/ 

♦define 

TO RAD 

1.74532925199e-2 

/♦Convert degrees to radians 

*/ 

♦define 

DEG30 

0.523598775598 

/* 

30 degrees in radians 


♦define 

DEG60 

1.04719755120 

/* 

60 degrees in radians 

*/ 

♦define 

DEG90 

1.57079632679 

/* 

90 degrees in radians 

*/ 

♦define 

DEG120 

2.09439510239 

/* 

120 degrees in radians 

*/ 

♦define 

DEG180 

3.14159265359 

/* 

180 degrees in radians 

♦/ 

♦define 

DEG270 

4.71238898038 

!* 

270 degrees in radians 

*/ 

♦define 

DEG360 

6.28318530718 

1 * 

360 degrees in radians 

*/ 


♦define 

J2 

1.08264e-3 

/* 

perturbation 

constant 

*/ 

♦define 

SUN_INCL 

-23.454*TO RAD 

/* 

pseudo-incl 

for plots 

*/ 

♦define 

SUN_TU 

83711.26224 

/* 

minutes per 

Solar TU 

*/ 


♦define CLASSICAL 1 /* parameter display options */ 

♦define SECONDARY 2 
♦define POSITIONS 3 
♦define OPTIONS 4 
♦define COMBINED 5 

♦define LETTER 18 /* font size and line spacing 

♦define LINE 22 

♦define ERROR -9999 

♦define max(x, y) (x > y) ? x : y /* return maximum value */ 

♦define nonzerop(x) fabs(x) > 0.00001 /* ~ float !“ 0.0 */ 

♦define zerop(x) fabs(x) < 0.00001 /* ~ float ==0.0 *! 
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•SMt.h Mon S«p 3 17:59:00 1990 3 

typedef struct { 

float xtriNEPOINTS]; 
float ylFINEPOINTSJ; 
float z[FINEFOINTS]; 
int pentFINEPOINTS]; 

) fine_clata; 


typedef struct ( 

float x[COARSEPOINTS]; 
float yCCOARSEPOINTS]; 
float z[COARSEPOIHTS]; 
int pen[COARSEPOIMTS]; 
) coarse data; 


typedef struct { 

float x[BORDERPOINTS]; 
float yfBOROERPOINTS) ; 
float 2 [BORDERPOINTS]; 
int pen[BORDERPOINTS]; 
) border_data; 


typedef struct ( 

char *name; 
int type; 
char *value; 
I Par; 


/* Displayable parameters */ 

/* for dialog boxes */ 

/* DISTANCE, SPEED, TIME */ 
/* cast this to whatever *! 


typedef struct { 

int numparams; 
Par *param3; 

) ParBlock; 


/* Block of parameters */ 

/* number of parameters *! 
/* list of parameters */ 


typedef struct { /* 

int orgx,orgy, 
sizex,sizey, 
gid; 

f window_struct; 


Window 

structure 

V 

/* 

lower left corner 

*/ 

/* 

window size 

*/ 

/* 

window grahics ID 

*/ 


typedef 


struct 

launch_st ruct 

{ /* 

Launch structure 

*/ 

char 

name f 3 0]; 

/* 

satellite name 

*/ 

float 

time. 

/* 

when burnout will occur 

*/ 


V, 

/* 

burnout velocity 



alt. 

/* 

burnout altitude 

*/ 


lat, 

1 * 

earth-fixed latitude 

*/ 


Ion, 

/* 

earth-fixed longitude 

*/ 


elev. 

/* 

elev angle - above horizon 

*/ 


azim; 

/* 

azimuth angle 

*/ 

struct 

launch struct 

*next 

;/* the next launch 

*/ 


) launch_struct ; 
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Mon Smp 3 17:59:00 1990 


typedef struct plot_struct ( /* Plot structure */ 
float lat,lon; /* a plot 

struct plot_struct *next; /* the next plot 
) plot_struct; 


struct 

roan_struct { /* 

Maneuver structure 

*/ 


int 

normal_maneuver 

;/* 

flag for deltav,yaw,pitch 

*/ 


float 

time. 

/* 

when maneuver will occur 

*/ 



anom. 

/* 

true anomaly of occurance 

*/ 



deltav. 

/* 

Change in velocity 

*/ 



yaw. 

/* 

yaw or heading change 

*/ 



pitch; 

/* 

pitch or flight path change 

*/ 


struct 

man_struct *next 

;/* 

the next maneuver 

*/ 


) man_ 

struct; 





struct 

satellite ( 

/* 

Satellite structure 


*/ 

char 

name[30]; 

/* 

satellite name 


*/ 

float 

major. 

/* 

semi-major axis of orbital 

ellipse 

*/ 


minor. 

/* 

semi-minor axis of orbital 

ellipse 

*/ 


ecc. 

/* 

eccentricity of orbital ellipse 

*/ 


incl. 

/* 

inclination of orbital ellipse 

*/ 


asc. 

/* 

longitude of the ascending 

node 

*/ 


peri. 

/* 

argument of periapsis 


*/ 


anom. 

/* 

true anomaly at epoch 


*/ 


epoch. 

/* 

time of epoch +ve (if after 

OOOOZ) 

*/ 


E, 

/* 

eccentric anomaly at epoch 


* / 


Pf 

/* 

X-coordinate in orbital plane 

*/ 



/* 

Y-coordinate in orbital plane 

*/ 


X, 

/* 

X-coordinate in IJK system 


*/ 


y> 

/* 

Y-coordinate in IJK system 


*/ 


z, 

/* 

Z-coordinate in IJK system 


*/ 


range. 

/■* 

Sat distance from center of 

earth 

*/ 


apogee range. 

/* 

Sat distance at apogee 


*/ 


perigee range. 

/* 

Sat distance at perigee 


*/ 


lat, 

/* 

Satellite latitude 


*/ 


Ion, 

/* 

Satellite longitude 


*/ 


period. 

/* 

Satellite period (one-revolution) 

*/ 


velocity. 

/* 

current satellite velocity 


*/ 


asc change. 

/* 

pertrubation constant 


*/ 


peri_change. 

/* 

pertrubation constant 


*/ 


user_floatl. 

/* 

for program modifications 


*/ 


user_float2. 

/* 

for program modifications 


*/ 


matrix[4][4]; 

/* 

orbit's rotation matrix 


*/ 

int 

plot. 

/* 

are we plotting ? 


*/ 


fill. 

/* 

do we fill in the orbit ellipse? 

*/ 


vectors. 

/* 

do we draw vectors on the e 

llipse? 

*/ 


color. 

/* 

display color 


*/ 


user inti. 

/* 

for program modifications 


*/ 


user int2; 

/* 

for program modifications 


*/ 


plot_struct *plots; /* 
man_struct *maneuver; /* 
struct satellite *next;/* 
) satellite; 


plots since plotting started 

next maneuver 

pointer to next satellite 
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SOMt .h 


Mon Sep 3 17:59:00 1990 


5 


/* global variable 

declarations 

*/ 


extern 

satellite header sat. 

/* 

Header to sat list 

*/ 

*current_sat 

; !* 

Current satellite 

*/ 

extern 

launch_struct *pending 

_launches;/* 

List of upcoming launches 

*/ 

extern 

float clock. 

/* 

Time elapsed 

*/ 

clock speed. 

/* 

time per minutes 

*/ 

viewer__x. 

/* 

Viewer X pos in Omniview 

*/ 

viewer_y. 

1* 

Viewer Y pos in Qmniview 

*/ 

viewer_z. 

/* 

Viewer Z pos in Omniview 

*/ 

viewer_range; 

/* 

Viewer Range in (Xnniview 

*/ 

extern 

int top menu. 

/* 

Menu ID 

*/ 

add_menu. 

/* 

Menu ID 

*/ 

orbit menu. 

/* 

Menu ID 

*/ 

color menu. 

/* 

Menu ID 

*/ 

picture menu. 

/* 

Menu ID 

*/ 

param menu. 

/* 

Menu ID 

*/ 

window order. 

/* 

which-where ex. OMNISATMERC 

*/ 

omni view gid. 

/* 

gid of OMNI view window 

*/ 

sat_view_gid. 

/* 

gid of SAT view window 

*/ 

merc_view_gid. 

/* 

gid of MERC view window 

*/ 

unit_type. 

!* 

CANONICAL, METRIC, AMERICAN 

*! 

param option. 

/* 

Which satellite parameters? 

*/ 

clock_running. 

/* 

is the clock running? 

*/ 

fineflag. 

/* 

Detailed Earth drawing? 

*/ 

borderllag. 

/* 

Draw political borders? 

*/ 

auto viewer. 

/* 

Auto pos for Omni view? 

*/ 

perturbations; 

/* 

Do we do perturbations? 

*/ 


/* Globe and map data records */ 

extern fine_data fine2d_data, fine3d_data; 
extern coarse_data coarse2d__data, coarse3d_data; 

extern border_data border2d_data, border3d_data; 

/* parameter block containing all displayable parauneters */ 

extern 

FarBlock classicalparamblock, 
secondaryparamblock, 
positionparamblock, 
optionsparamblock, 
combinedparamblock, 
sysparaunblock; 
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smet .h 


Mon S«p 3 17:59:00 1990 


6 


extern 

char *instructions[28], 

*color_naine3 [MAX_COLORS] 
*unit_nan\e3 [3], 
•speed_unit[3], 
*clistance_unit [3]; 

extern 

float speed_factor[3], 

di3tance_factor[3]; 

extern /* 

window_struct backwin, 
mainwin, 
auxwinl, 
auxwin2, 

3atli3twin, 

dialogwin, 

paramwin; 


/* 

User 

instructions 

*/ 

, /* 

Nai.ies 

of 

all 

the colors 

*/ 

/* 

Names 

of 

the 

Units 

*/ 

/* 

Names 

of 

the 

speed units 

*/ 

/* 

Names 

of 

the 

distance units 

*/ 


/* Speed conversion factor 

*/ 

/* Distance 

conversion factor 

*/ 

up all our 

window parameters 

*/ 


extern 

float ident_matrix[4] [4] , 

Color[MAX COLORS] (4] ; 


/* Nice constant to have laying around */ 
/* Array of colors for c3f or c4f */ 


/* forward function declarations */ 

draw_instructions(); 
add_launch(); 
ask_and_read_sat_file(); 
a3k_and_write_sat_file(); 
add_maneuver(); 
add_3ateliite(int); 
cu3tom_3atellite() ; 
delete_3atellite() ; 
take_picture() ; 
program_exit() ; 
smet_param_option(int); 

3met_3at_color(int); 

float ask time(); 
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It vars. c 


Mon Sap 3 19:00:52 1990 


1 


/***************^**»************** 

/* 

/* 9met_var3.c 

/* by Carlos I. Noriega 

/* 

/* This file contains the global variables, constants and 

/* structures used by the SMET demonstrator. 

/* 

/* Written - September 3, 1990 

/* Modified - September 3, 1990 

/* 

/********************************* 


*! 

*/ 

*/ 

*! 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 


♦include "smet.h" 


satellite header sat. 

/* 

Header to sat list 

*/ 


*current_sat — NULL; 

/* 

Pointer to list of sats 

*/ 

launch struct Spending launches 

NULL;/* Upcoming launches 

*/ 

float 

clock ” 0.0, 

/* 

Time elapsed in minutes 

*/ 


clock_3peed - 30.0, 

/* 

Time per cycle in minutes 

*/ 


viewer x. 

/* 

Viewer X pos in Omniview 

*/ 


viewer y. 

/* 

Viewer T pos in Cwnniview 

*/ 


viewer_z. 

/* 

Viewer Z pos in Omniview 

*/ 


viewer range; 

/* 

Viewer Range in Omniview 

*/ 

int 

top menu. 

/* 

Menu ID 

*/ 


add_menu. 

/* 

Menu ID 

*/ 


orbit_menu. 

/* 

Menu ID 

*/ 


picture_menu. 

/* 

Menu ID 

*/ 


color menu. 

/* 

Menu ID 

*/ 


param menu. 

/* 

Menu ID 

*/ 


window order. 

/* 

which-where ex. C^INISATMERC 

*/ 


omni view gid. 

/* 

gid of OMNI view window 

*/ 


sat view gid. 

/* 

gid of SAT view window 

*/ 


mere view gid. 

/* 

gid of MERC view window 

*/ 


unit type ” CANONICAL, 

/* 

CANONICAL, METRIC, ENGLISH 

*/ 


param option “ CLASSICAL, 

/* 

Which satellite parameters? 

*/ 


clock running = TRUE, 

/• 

is the clock running? 

*/ 


fineflag = TRUE, 

/* 

Detailed Earth drawing? 

*/ 


borderflag - TRUE, 

/* 

Draw political borders? 

*/ 


auto viewer *• TRUE, 

/* 

Auto pos for Omni view? 

*/ 


perturbations “ EALSE; 

/* 

Do we do perturbations? 

*/ 


/* Globe and map data records 

fine_data fine2d_data, fine3d_data; 
coar3C_data coatse2d_data, coarse3d_data; 
border data border2d data, border3d_data; 
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siMt_Tars.c Mon Sap 3 19:00:52 1990 

/* data structures for parameter display */ 


2 


Tar classical_list[] = ( 

( "Name ", STRING, (char *)header_sat.name ), 

( "Semi-Major ", DISTANCE, (char *)&(header_sat.major) ), 

( "Semi-Minor ", DISTANCE, (char ")S(header_sat.minor) ), 

( "Eccentricity ", FLOAT, (char *)4(header_sat.ecc) ), 

( "Inclination ", DEGREE, (char *)fi(header_sat.incl) }, 

{ "Ascending Node ", DEGREE, (char ")4(header_sat.asc) ), 

( "Arg of Periapsis ", DEGREE, (char *)4(header_sat.peri) ), 

( "True Anomaly ", DEGREE, (char *)4(header_sat.anom) ), 

{ "Epoch ", TIME, (char *)4(header_aat.epoch) ) ), 

secondary_list{] = { 

{ "Name ", STRING, (char *)header_sat.name ), 

( "Perigee ", DISTANCE, (char ")4(header_sat.perigee_range) ), 
{ "Apogee ", DISTANCE, (char *)4(header_sat.apogee_range) ), 

( "Epoch ", TIME, (char ")4(header_3at.epoch) ), 

{ "Period ", TIME, (char *)4(header_sat.period) ), 

( "Range ", DISTANCE, (char *)4(header_3at.range) ], 

{ "Velocity ", SPEED, (char *)4(headcr_3at.velocity) ) ), 


po3ition_li3t[1 = ( 


( 

"Name 

", STRING, 

(char ")header 

sat.name 

) 

{ 

"X 

n 

f 

DISTANCE, 

(char *)4(header 

sat .x) 

K 

{ 

"Y 

n 

t 

DISTAIJCE, 

(char *)4(header 

_sat.y) 

), 

{ 

"Z 

ft 

/ 

DISTANCE, 

(char *)4(header 

_sat.z) 

), 

{ 

"P 

1* 

f 

DISTANCE, 

(char *)4(header 

_sat.p) 


{ 

"Q 

« 

! 

DISTANCE, 

(char *)4(header 

_3at.q) 

)r 


{ "Latitude ", LAT, (char *)4(header_sat.lat) ), 

( "Longitude ", LON, (char *)4(header_3at.Ion) ) ), 


option_list[) ■ { 

( "Name ", STRING, (char ")header_sat.name }, 

( "Color ", COLORCODE, (char *)4(header_sat.color) ), 
{ "Fill ", FLAG, (char *)4(header_sat.fill) ), 

( "Vectors ", FLAG, (char *)4(header_sat.vectors) ), 

( "Plot ", FLAG, (char *)4(header_sat.plot) ) }, 


sysparamlist[] ” ( 


{ 

( 

I 

( 

{ 

( 

{ 

( 

( 

( 

( 


"Clock ", FLAG, (char *)4clock_running ), 

"Current Time ", TIME, (char *)4clock }, 

"Clock Speed ", TIME, (char *)4clock_3peed ), 
"Detailed Maps ", FLAG, (char *)4fineflag ), 
"Political Borders ", FLAG, (char *)4borderflag ), 
"Perturbations ", FLAG, (char *)4petturbations ), 
"Auto Omniviewer ", FLAG, (char *)4auto_viewer ), 
"Viewer X ", DISTANCE, (char ")4viewer_x ), 

" Y ", DISTANCE, (char *)4viewer_y ), 

" Z ", DISTANCE, (char *)4viewer_z ), 

"Units ", UNITCODE, (char ")4unit_type ), 
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it Tatra, c 


Mon S«p 3 19:00:52 1990 


3 


combined_list n ” ( 

( "Name ", STRING, (char *)header_sat.name ), 

( "Semi-Major ", DISTANCE, (char *)fi(header_sat.major) ), 

( "Semi-Minor ", DISTANCE, (char ")&(header_sat.minor) }, 

{ "Eccentricity ", FLOAT, (char *)S{header_3at.ecc) ), 

{ "Inclination ", DEGREE, (char ")i(header_sat.incl) ), 

( "Ascending Node ", DEGREE, (char ")&(header_sat.asc) ), 

{ "Arg of Periapsis ", DEGREE, (char *)&(header_3at.peri) ), 

( "True Anomaly ", DEGREE, (char *)£(header_3at.anom) ), 

( "Epoch ", TIME, (char *)£(header_sat.epoch) ), 

{ "Period ", TIME, (char ")£(header_3at.period) ), 

{ "Apogee ", DISTANCE, (char *)£(header_3at.apogee_range) }, 

( "Perigee ", DISTANCE, (char *)£(header_sat.perigee_range) ), 
{ "Range ", DISTANCE, (char *)£(header_sat.range) ), 

( "Velocity ", SPEED, (char *)£(header_3at.velocity) }, 


( 

"X 

fl 

r 

DISTANCE, 

(char 

*)£(header_ 

sat.x) 

) 

{ 

"Y 

m 

f 

DISTANCE, 

(char 

*)£(header_ 

sat. y) 

) 

( 

"Z 

fl 

f 

DISTANCE, 

(char 

")£(header_ 

sat.z) 

) 

{ 

"P 

n 

t 

DISTANCE, 

(char 

")£(header_ 

sat.p) 

) 

1 

"Q 

If 

t 

DISTANCE, 

(char 

*)£ (header_ 

sat .q) 

) 


( "Latitude ", LAT, (char *)fi(header_3at.lat) ), 

{ "Longitude ", LON, (char *)fi(header_3at.Ion) ), 

( "Color ", COLORCODE, (char *)£(header_sat.color) ), 
( "Fill ", FLAG, (char *)£(header_sat.fill) ), 

( "Vectors ", FLAG, (char *)£(header_3at-vectors) ), 

{ "Plot ", FLAG, (char *)£(header_sat.plot) ) ) ); 


ParBlock classicalparamblocic « { 
9, 

classical_list 

), 

secondaryparamblock »■ { 

7, 

secondary_list 

), 

positionparamblock — ( 

8 , 

position_liSt 

), 

optionsparamblock « ( 

5, 

option_list 

», 

corobinedparamblock — ( 

25, 

combined_list 

), 

sysparamblock - { 

11 , 

sysparamlist 

}; 


/* 

ntanparams 

*/ 

/* 

Par list 

*/ 

!* 

numparams 

*/ 

/* 

Par list 

*/ 

/* 

numparams 

*/ 

/* 

Par list 


/* 

numparams 

*/ 

/* 

Par list 

*/ 

/* 

numparams 

*/ 

/* 

Par list 

V 

/* 

numparams 

*/ 

/* 

Par list 

*/ 


62 




it vars.c 


Mon S«p 3 19:00:52 1990 


4 


char *instruction3[28] - /* User instructions */ 

{" SATELLITE MANEUVER EVALUATION TOOL" 

"This program is designed to let the user view how the manipulation", 
"of different orbital parameters affect the path of a satellite.", 

"You will see the orbiting satellite from an omniscient point of", 
"point of view, the earth from the satellite's point of view and a", 
"mercator projection map with the track of the satellite marked.", 

"The user can also launch or maneuver satellites by use of a menu", 
"command.", 

If 19 91 19 

f f 

"INSTRUCTIONS 

"Rightmouse: Used to pop up menus to execute the desired command.", 

99 99 

t 

19 19 

f 

"Middlemouse: Used to select a parameter, a satellite or an ", 

" Auxiliary Window.", 

99 99 

r 

"Dials: The six upper dials will allow the user to adjust following", 

" parameters. The order is: SEMI-MAJOR AXIS ECCENTRICITY", 

" INCLINATION ASCENDING NODE", 

" PERIGEE CLOCK SPEED", 

n n n 99 
r t 

"Press the Rightmouse button to continue."), 

*colorenames[MAX_COLORS] - { "Black", "Red", "Green", "Yellow", "Blue", 

"Magenta", "Cyan", "White","Earth", "Coast", 
"Border", "LatLong", "Sun","Red", "Green", 
"Yellow", "Blue", "Magenta", "Cyan", "White", 
"Gray"), 

*unit_name3[3) = { "Canonical", "Metric", "English" }, 

*3peec'_unit [3] - { "DU/TU", ")an/sec", "ft/sec" ), 

*distance_unit[3] - { "DU", "km", "NM"); 

float speed_factor[3] - { 1.0, DU_PER_TU_METRIC, DU_PER_TU_ENGLISH ), 
distance_factor[3] - { 1.0, DU_METRIC, DU ENGLISH ); 


63 






•Mt ▼ars.c Mon S«p 3 19:00:52 1990 5 

/* Set up all our window parameters allowing for transfer */ 

/* to a different size screen system */ 

window struct /* lower-left corner, window-size and dummy GID •/ 


backwin = { 0, 0, XMAXSCREEN, YMAXSCREEN, 0 ), 


mainwin 

- ( 

XMAXSCREEN 

* 

265/1279, 

XMAXSCREEN 

* 

273/1023, 





XMAXSCREEN 

it 

750/1279, 

XMAXSCREEN 

* 

750/1279, 

0 

), 

auxwinl 

- ( 

XMAXSCREEN 

it 

315/1279, 

XMAXSCREEN 

* 

23/1023, 





XMAXSCREEN 

it 

250/1279, 

XMAXSCREEN 

it 

250/1279, 

0 

), 

auxwin2 

- ( 

XMAXSCREEN 

it 

715/1279, 

XMAXSCREEN 

it 

23/1023, 





XMAXSCREEN 

it 

250/1279, 

XMAXSCREEN 

* 

250/1279, 

0 

), 


satlistwin - { XMAXSCREEN * 1015/1279, YMAXSCREEM * 0/1023, 

XMAXSCREEN * 264/1279, XMAXSCREEN * 1023/1023, 0 ), 

dialogwin - ( XMAXSCREEN * 270/1279, XMAXSCREEN * 620/1023, 

XMAXSCREEN * 740/1279, XMAXSCREEN * 45/1023, 0 ), 

paramwin = ( XMAXSCREEN * 0/1279, XMAXSCREEN * 0/1023, 

XMAXSCREEN * 265/1279, XMAXSCREEN * 1023/1023, 0 ); 


float ident matrix(4](4] - ( (1.0, 

0.0, 

0.0, 

0.0) 

f 

(0.0, 

1.0, 

0.0, 

0.0) 

f 

(0.0, 

0.0, 

1.0, 

0.0) 

f 

(0.0, 

0.0, 

0.0, 

1.0) 

) 


Color (MAX_COLOF.3] f 4] - {(0.000, 0.000, 0.000, 0.000), 

(1.000, 0.000, 0.000, 0.000), 
(0.000, 1.000, 0.000, 0.000), 
(1.000, 1.000, 0.000, 0.000), 
(0.000, 0.000, 1.000, 0.000), 
(1.000, 0.000, 1.000, 0.000), 
(0.000, 1.000, 1.000, 0.000), 
(1.000, 1.000, 1.000, 0.000), 
(0.250, 0.800, 0.500, 0.000), 
(0.750, 0.200, 0.000, 0.000), 
(0.500, 0.250, 0.000, 0.000), 
(0.750, 0.750, 0.750, 0.000), 
(1.000, 0.800, 0.200, 0.400), 
(0.750, 0.000, 0.000, 0.200), 
(0.000, 1.000, 0.000, 0.200), 
(1.000, 1.000, 0.000, 0.200), 
(0.000, 0.500, 1.000, 0.200), 
(1.000, 0.000, 1.000, 0.200), 
(0.000, 1.000, ’.000, 0.200), 
(1.000, 1.000, 1.000, 0.200), 
(0.600, 0.600, 0.600, 0.200) ); 


/* 

BLACK 

*/ 

/* 

RED 

•/ 

/* 

GREEN 

*/ 

/* 

XELLOW 

*/ 

/* 

BLUE 

*/ 

/* 

MAGENTA 

*! 

/* 

CXAN 

*/ 

/* 

WHITE 

*/ 

/* 

EARTH 

*/ 

/* 

COAST 

*/ 

/* 

BORDER 

*/ 

/* 

LATLONG 

*/ 

/* 

SUN 

*/ 

/* 

SAT_RED 

*/ 

/* 

SAT_GREEN 

*/ 

/* 

SAT_XELLOW 

*/ 

/* 

SAT BLUE 

*/ 

/* 

SAT MAGENTA 

*/ 

/* 

SAT CXAN 

*/ 

/* 

SAT_WHITE 

*/ 

/* 

SAT GRAX 

*/ 
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y* *^k*'************^******^**'****t*** ★/ 

/* */ 

/* 3met_addelete.c */ 

/* by Carlos I. Noriega */ 

/* */ 

/* This file contains the routines used by SMET to */ 

/* add and delete satellites, add and delete plots and */ 

/* switch and crash satellites. */ 

/* */ 

/* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/* */ 

/A*********************************/ 

iinclude "smet.h" 

/* MAKE_NEW_SATELLITE */ 

/* Creates a new satellite and sets all default settings *! 

/* Asks user for name if none is provided. */ 

make_new_current_3at(name) 
char *name; 

( 

satellite *3at; 

sat = (satellite *)malloc(sizeof(satellite)); 

3at->epoch - clock; 

if (strlen(name)) strcpy(sat->name,name); 

else dialog ("Enter name of new satellite.", sat->name); 

3at->color - SAT_RED + rand() % 6; 

3at->fill - FALSE; 

3at->vector3 - FALSE; 
sat->plot “ FALSE; 
sat->plot3 •= NULL; 
sat->maneuver » NULL; 

3at->next •- current_3at; 
header_sat next “ current_3at — sat; 
u3er_new_sat(); 

) 
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/* ADD_SATELLITE */ 

I* Set the five basic parameters of the orbit and init- */ 

/* ialize the dials as appropriate. */ 

add_satellite(n) 

int n; /* orbit type selected by pop up menu •/ 

( 

make_ne»*_current_sat ("") ; 
switch (n)( 


case 1: /* 300 nm circular LEO */ 

current_aat->major — 1 + 300/distance_factorlENGLISH]; 
current_sat->ecc • 0.0; 
current_sat->incl •» DEG30; 
current_sat->asc “ DEG30; 
current_sat->peri ■ 0.0; 
break; 

case 2: /* Molniya Orbit */ 

cutrent_sat->major — pow(ROTATION_RATE*TU*2, -2.0/3.0); 
current_sat->ecc “ 0.75; 
current_sat->incl — asin(sqrt(0.8)); 
current_sat->a3c - DEG60; 
cui.rent_sat->peri “ DEG270; 
break; 

case 3: /* Geosynchronous Orbit */ 

current_sat->maior = pow (ROTATION_RATE*TU, -2.0/3.0); 
current_3at->ecc “ 0.0; 
cuttent_sat->incl ” 0.0; 
current_3at->a3C * 0.0; 
current_3at->peri - 0.0; 
break; 

case 4: /* Custom Orbit */ 

ask_ma jor(); 
ask_ecc(); 
a3k_incl<) ; 
a3k_asc() ; 
a3k_peri() ; 
break; 

) 

update_static_paranieter3 (current_sat) ; 

update_sat_matrix(current_sat); 

reset dials (); /* set dials to current satellite 


*/ 
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/* DELETE_SATELLITE 

/* Removes the current satellite from SMET 
clelete_satellite () 

I 

header_aat.next - current_sat->next; 
delete_plot3(current_sat); 
free((char *) current_sat); 
current_sat - header_sat.next; 
if (current sat !“ MULL) reset_dials(); 


*/ 

*/ 


/* ADD_PLOT & DELETE_PLOTS */ 

/* Add a plot to the linked list or purge the list. */ 

add_plot(sat) 
satellite *sat; 

( 

plot_struct *new_plot; /* link for new lat/lon */ 

new_plot - (plot_struct *)malloc(sizeof(plot_3truct)); 
new_plot->lat ” 3at->lat; 
new_plot->lon «• sat->lon; 
rew_plot->next “ 3at->plot3; 

3at->plots - new_plot; 


delete plots (sat) 
satellite *sat; 

( 

plot_struct *temp; 

for (temp -> sat->plots; sat->plots !- NULL; temp - sat->plots) 
( 

sat->plot3 " sat->plot3->next; 
free((char *) temp); 

) 

) 
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/* SWITCH_CURRENT_SAT */ 

/* Figures out which satellite was chosen from mouse y */ 

J* position. Walk down list until you got it and switch'em. */ 

switch_current_sat(y) 
short y; 

int sat_number " (float)(satlistwin.sizey - y) / (float) LINE - 3, 
count; /* duimiy loop counter */ 
satellite *walker; 

if (aat_number > 1) 

( 

walker — current_sat; 

for (count —I; walker !— NULL tt count + 1 < s_'_number; count++) 
walker ■■ walker->next; 
set_next_to_current(walker) ; 

) 

) 


/* CRASH_SATELLITE */ 

/* Set the desired sat to the current position so that */ 

/* you can do a delete satellite. Called when range >-1. */ 

cra3h_3atellite(3at, previou3_3at) 
satellite *sat, ‘previous sat; 

( 

chat att(501; 

sprintf(str, "%s crashed into the earth and was deleted.",sat->name); 
message(str); 

set_tiext_to_cutrent (previous_3at) ; 
delete_satellite(); 
sat - previou3_3at; 


/* SET_NEXT_TO_CURRENT */ 

/* Used by two previous procedures to switch the linked */ 

/* list of satellites. */ 

3et_next_to_current(previous) 
satellite "previous; 

( 

if (previous !« NULL && previou3->next !“ NULL previous !■« (header_sat) 
( 

header_sat.next - previou3->next; 
previous->next » pteviou3->next->next; 
header_S8t.next-’next - current_sat; 
curtent_sat - header_sat.next; 
reset_dials(); 

) 

) 
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/*******************•*********«****/ 
/* */ 

/* smet_axe9.c */ 

/* by Carlos 1. Noriega */ 

/* ’*/ 

/* This file contains the routines to draw the axes */ 

/* used by the SMET demonstrator. */ 

/* */ 

/* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/* */ 

y^4r*A*A*AAA*A*A***A*A*****A******A**/ 


♦include "smet.h" 


/* DRAW-AXES */ 
/* draw the IJK axes at the center of the */ 
/* world coordinate system. Remember that */ 
/* IJK is not aligned the same as the systems */ 
/* XYZ coordinate system. */ 


draw_axcs() 

( 

/* go to white */ 
c3f (Color (WHITE]) ; 

/* set the linewidth - what a comment !! */ 

linewidth(3); 

/* label the I axis */ 
cmov(0.0, 0.2, 1.2); 
fmprstr("I"); 

/* label the J axis */ 
cmov(1.2, 0.2, 0.0); 
fmprstr("J" ) ; 

/* label the K axis */ 
cmov(0.2, 1.2, 0,0); 

fmprstr( "K") ; 

/* draw the J axis */ 
move (1.0, 0.0, 0.0); 
draw (12.5, 0.0, 0.0); 

/* draw the K axis */ 
move (0,0, 1.0, 0.0); 

draw(0.0, 12.5, 0.0); 

/* draw the I axis */ 
move (0.0, 0.0, 1.0); 
draw(0.0, 0.0, 12.5); 
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/*****************************★****/ 
/* */ 

/* smet_dialog.c */ 

/* by Carlos I. Noriega */ 

/* */ 

/* This file contains the routines to give the user messages */ 

I* and ask him questions. It is used by the SMET */ 

/* demonstrator. */ 

/* */ 

/* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/* */ 

/***************************«***** 

iinclude "smet.h" 

/* DIALOG & MESSAGE */ 

I* Use the dialogbo:: routines to send a message or ask a *l 

/* question. */ 


dialog(prompt,response) 
char ^prompt, *respon3e; 

( 

open_dialogbox(); 

get_user_response(prompt,response); 
clo3e_diaiogbox(); 

I 


message(prompt) 
char *prompt; 

\ 

open_dialogbox0; 
clearprompt(prompt); 
ringbell(); 
sleep (3); 

clo3e_dialogbox (); 
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DIALOGBOX Routines */ 

/* Uses the hidden dialogwin for all user messages. */ 

open_dialogbox() /* "POP" the box. */ 

( 

winset(dialogwin.gid); 

octho2(0.0,dialogwin.sizex, 0.0, dialogwin.sizey); 

loadmatrix(ident_matrix); 

singlebuffer(); 

gconfig() ; 

winpop(); 

zciear () ; 

qdevice(KEYED); 


clearprompt(prompt) /* Clear the box and type the message. */ 

char "prompt; 

( 

c3f(Color[CYAN]); 
clear(); 

c3f(Color[BLACK]); 
linewidth(2) ; 

recti(7, 7, dialogwin.sizex-7, dialogwin.si2ey-7); 

cn>ov2i(12, 20); 

charstr(prompt); char3tr(" "); 

» 


clo3e_dialogbox() /* "PUSH" the box. */ 

t 

zciear(); 
winpuslj () ; 
unqdevice(KEYED); 
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/* GET_USER_RESPONSE */ 

/* A simple parser for user interaction. Does not currently */ 

/* check if response string is shorter than fixed length of 40.*/ 

/* Ignores events other than the keyboard. */ 

get_user_response(prompt,response) 
char *prompt, *response; 

( 

int cur_str_len, 

donetyping - FALSE; 
short c; 
char *str; 

Device dev; 

responsetOJ — '\0'; 
cur_str_len — strlen(response) ; 

clearprompt(prompt); /* display prompt */ 

while(!donetyping) /* read till eof ('\n' or '\r') */ 

if(dev — qread(&c)) 

( 

if(dev !- KEYED) continue; /* don't care */ 

switch(c) 

( 

case '\027': /* sets cursor back to start */ 

responsetO] •• ' \0'; 
cur_3tr_len ■■ strlen (response) ; 
clearprompt(prompt); 
break : 

case ' \n': 

case '\r': /* carriage return or enter */ 

donetyping — TRUE; 
break; 

case 127: 

case '\b': /* backspace or delete */ 

if(cur_stt_len) ( 

response!—cur_str_lenj — '\0'; 
clearprompt(prompt); 
charstr(response); ) 

else ringbellO; 
break; 

default: /* add to string */ 

if(cur_str_len < 41) ( 

str - 4response(cur 3tr_len]; 
response[cur_str_len++] — c; 
response[cur_str_len] — '\0'; 
charstr(str); ) 

else ringbellO; 
break; 

) 

I 

I 
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/* ******************************** *^ 
/* */ 

/* smet_files-c */ 

/* by Carlos I. Horiega •/ 

/* */ 

/* This file contains the routines used by SMET to */ 

/* read and write files. *! 

/* */ 

/* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/* */ 

/**********************************^ 


♦include "smet.h" 


/* READ_WORLD_FILE 

*/ 

/* Reads earth data for map and globe from file in the 

*/ 

/* current directory. 


read world fileO 


\ 

FILE *fd; 


fd * fopen("3met data.d","r"); 


if(fd -= NULL) 



( 

message("Unable to open file 3met_data.d, quiting SMET"); 
program_exit(); 

) 

else 

« 

Iseek(fd,OL,0); 

fread(ifine3d_data, sizeof(fineSd data), 1, fd); 
fread(ifine2d_data, sizeof(fine2d_data), 1, fd); 
fread(icoarse3d_data, sizeof(coarse3d_data), 1, fd); 
fread(icoar3e2d_data, sizeof(coar3e2d_data), 1, fd); 
fread(iborder3d_data, sizeof(border3d_data), 1, fd) ; 
fread(iborder2d_data, sizeof(bordor2d_data), 1, fd); 
fclose(fd); 

) 

) 
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/* ASK_AND_READ_SAT_FILE & READ_SAT_FILE */ 

/* Reads satellite data previously saved. */ 

a3k_and_read_sat_file() 

( 

char name[40]; 

dialog("Enter the name of the satellite filename); 
read sat file(name); 


read_3at_file(file_name) 
char file_name(]; 

{ 

FILE *fd; 
satellite "sat; 
int magic; 
char str[50); 

fd = fopen(file_name,"r"); 
if(fd -- NULL) 

{ 

sprintf(str,"Unable to open file - %3.",file_name); 
message(str); 

) 

else 

{ 

Iseek(fd,OL,0); 

fread (&n\agic, sizeof(int), 1, fd) ; 

if (magic !- 4567) /* ensure properly formatted file */ 

{ 

sprintf(str, "%s is not a satellite file.file_name); 
message(str); 

) 

else 

{ /* read satel)rtes until you get a read error */ 

sat = (satelli’-e *)malloc (sizeof (satellite) ) ; 
while (fread(3at, sizeof (satellite), 1, fd)) 

I 

3at->plots “ NULL; 

3at->maneuver — NULL; 
sat->next •• current_sat; 
current_sat - sat; 

sat — (satellite *)malloc(sizeof(satellite)); 

) 

free((char *) sat); 

header_3at.next ■= current_sat; /* make last one current "/ 
reset_dial3(); 

I 

fclose(fd); 

) 

) 
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/* ASK_AND_WRITE_SAT_FILE & WRITE_SAT_FILE */ 

/* Write satellite data into a file. */ 

as)c_and_write_sat_f ile () 

{ 

char name[40], number[5]; 

dialog("Enter the name of the satellite file.", name); 
dialog("Enter the number of the satellites saved.", number); 
write sat file(name,atoi(number)); 


write_sat_f ile (file_naune, numsats) 
char file_name[] ; 
int numsats; 

( 

FILE *fd; 
satellite *3at; 
int c, magic = 4567; 
char str [50] ; 

fd = fopen (f ile_nam\e, "w") ; 
if(fd == NULL) 

( 

sprintf(str,"Unable to open file - %s.",file_name); 
message(str); 

) 

else 

( 

Iseek(fd,OL,0); 

fwrite(imagic, sizeof(int), 1, fd) ; /* ensure properly formatted file •/ 
sat “ current_3at; 

for (c •• 0; (c < numsats 4S sat • ” NULL); C++) 

( 

fwrite(sat, sizeof(satellite), 1, fd); 
sat “ 3at->next; 

) 

fclose(fd); 

} 

) 
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/* TAKE-PICTURE */ 
/* Use image library routines to write an image */ 
/* in the selected format. Black and White are */ 
/* inverte-^ *or use in EPS format. */ 


take_picture(num) 
int num; 

( 

NPSIMAGE *imgl,*img2, *img3; /* temp ptr to images */ 
char filename[50], 
str[751; 

switch (num) 

I 

case 1: 

dialog("Enter the name for the RGB file filename); 

message("Saving the screen will take about 90 seconds."); 
refresh_main(); 

save_window_as_sgi_rgbimage(filename,backwin.gid); 
break; 
case 2: 

dialog("Enter the name for the RGB file filename); 

message("Saving the screen will take about 90 seconds."); 
invert_bw(); 

save_window_a3_sgi_rgbimage(filename, backwin.gid); 
revert_bw () ; 
break; 
case 3: 

dialog("Enter the name for the EPS file -", filename); 
message("Saving the screen will take about 60 seconds."); 
invert_bw(); 

imgl = read_window_fb_as_rgb_npsimage(backwin.gid,filename) ; 
img2 » 3hrink_np3image(imgl, 2); 
delete_an_image(imgl); 
img3 - rotate_np3image(img2); 
delete_an_iniage (img2) ; 
write_eps_image(filename, img3, 125); 
delete_an_image(img3) ; 
revert_bw(); 
break; 
case 4: 

dialog("Enter the name for the EPS file filename); 
message("Saving the main window will take about 30 seconds.") 
invert_bw(); 

imgl - read_window_fb_as_rgb_npsimage(mainwin.gid, filename); 

img2 - shrink_npsimage(imgl, 2); 

deiete_an_image(imgl) ; 

write_eps_image(filename, img2, 125); 

delete_an_image(img2) ; 

revert_bw(); 

break; 

) 

) 
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*/ 
*/ 


/* Support routines for TAKE_PICTURE 


ref resh_nvain ') 

{ 

switch(window_order) /* refresh main window after message */ 

( 

case OMNISATMERC: 

case OMNIMERCSAT: draw_omni_view(); break; 
case SATCMNIMERC: 

case SATMERCOMNI: draw_sat_view(); break; 
case MERCOMNISAT; 

case MERCSATOMNI: draw_merc_view(); break; 

) 

) 


invert_bw() 

( 

Color IBLACK] [0] - Color[BLACK] [1] - Color[BLACK] [2] - Color[BLACK] [3] - 1.0; 
Color[WHITE][0] - Color[WHITE][1] - Color[WHITE][2] - Color[WHITE][3] - 
Color[SAT_WHITE][0] - Color[SAT_WHITE1[1] - Color[SAT_WHITE][2] - 

Color[SAT_WHITE][3] - 0.25 

winset(backwin.gid) ; 
c3f(Color[BLACK]); 
clear() ; 
swapbuffers() ; 
draw_window3() ; 

draw windows 0; /* draw twice to ensure that proper buffers are showing */ 


revert_bw() 

( 

Color[BLACK][0] - Color[BLACK][1] - Color[BLACK][2] - Color[BLACK][3] - 0.0; 
Color[WHITE][0] - Color[WHITE][1] - Color[WHITE][2] - Color[WHITE][3] - 
Color(SAT_WHITE][0] - Color[SAT_WHITE][1] - Color[SAT_WHITE)[2] - 

Color[SAT_WHITE)[3] - 1.0; 

winset(backwin.gid); 
c3f(Color[ black]); 
clear (); 
swapbuffers(); 
draw windows() ; 
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/**********************************/ 


/* */ 

/* sinet_globe. c */ 

/* by Carlos I. Noriega */ 

/* *l 

/* This file contains the routines used by SMET to */ 

/* draw the globe in a viewing window. */ 

/* */ 

/* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/’» */ 


******************************** */ 
♦include "smet.h" 

/* DRAW-GLOBE */ 

/* Assumes that a central disk normal to viewer drawn *! 

/* by calling routine. */ 

draw_globe(rotation) 

float rotation; /* amount of y-rotation of the globe */ 

[ 

pushmatrix(); 

rot(rotation*TO_DEG,'y'); /* turn the globe */ 

draw_coasts3d(); 
draw_borders3d0 ; 
draw_latlons3d() ; 

popmatrix() ; 

I 


/* DRAW-COASTS 3D */ 

/* Draw the desired data TINE or COURSE. */ 

draw_coasts3d () 

( 

int c; 

linewidth (2); 
c3f(Color[COAST]); 
if (fineflag) 

for (c - 0; c < FINEPOINTS; C++) 
if (fine3d_data.pen[c]) 

draw(fine3d_data.X[c], fine3d_data.y(c], fine3d_data.z[c]); 
else 

move(fine3d_data.x(c], fine3d_data.y[c], fine3d_data.z[c]); 

else 

for (c - 0; c < COARSEPOINTS; C++) 
if (coarse3d_data.pen[cj) 

draw(coarse3d__data.x(cl, coar3e3d_data.y[cj, coarse3d_data.z[c]) 
else 

move(coarse3d data.x(c], coarse3d_data.y(c], coarse3d_data.z[c]) 

» 
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*/ 

*/ 


/* DRAW-BORDERS3D 

/* Draw the border data if desired. 

draw_borders3d() 

{ 

int c; 

linewidth(1) ; 
c3f(Color[BORDER]); 
if (borderflag] 

for (c - 0; c < BORDERPOINTS; C++) 
if (border3d_data.pen[c]) 

draw(border3d_data.x[c], border3d_data.y[c], border3d_data.z[c]) 
else 

move(border3d data.x[c], border3d data.y[c], border3d data.z[c]) 


/•* DRAW-LATLONS3D */ 

/* Draw the latitude and longitude lines on the globe. */ 

draw_latlon33d() 

( 

int c; 

linewidth(1); 

c3f(Color[LATLONG]); 

pushmatrix(); 

/* draw the longitude lines */ 
pushmatrix(); 

for (c “ 1; c <* 12; C++) /* draw an arc every 30 degrees */ 

1 

rotate(300, 'y'); 

arc(0.0,0.0, 1.0,-900, 900) ; 

) 

popmatrix() ; 


/* draw the latitude lines */ 

rotate(900,'x'); /* flip over on its side */ 

circ (0.0,0.0,1.0); /* equator */ 

translate(0.0,0.0,-0.5); 

circ(0.0,0.0,0.866); /* 30 degrees north */ 

translate(0.0,0.0,-0.366); 

cite (0.0,0.0,0.5); /* 60 degrees north */ 

translate(0.0,0.0, 1.366) ; 

circ(0.0,0.0,0.866); /* 30 degrees south */ 

translate(0.0,0.0,0.366); 

circ(0.0,0.0,0.5); /* 60 degrees south */ 


popmatrix(); 

) 
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/********i»*************************/ 

/* *! 

/* sinet_axes.c */ 

/* by Carlos I. Noriega */ 

/* */ 

/* This file contains the initialization routines */ 

/* used by the SMET demonstrator. */ 

/* */ 

I* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/* */ 

/********************»*************/ 


finclude "smet.h" 


I* INITIALIZE_SMET */ 
/* Open windows, set up menus, initialize queue, */ 
/* draw the user instructions. */ 


initialize_smet(argc,argv) 
int argc; 
char **argv; 

( 

int c; 

fmfonthandle param_font; 

/* initialize windows 

init_window(fibackwin) ; 
init_window(idialogwin); 
init_window(tsatlistwin) ; 
init_window(&paramwin); 
init_window(Smainwin); 
init_window(sauxwinl); 
init_window (Savixwin2) ; 

message (" 

/* make the pop up menus 


*/ 


/* a black background */ 
/* the dialog window */ 
/* where the satellite names are listed */ 
I* where the parameters are listed */ 
/* the centered large window */ 
/* the lower left small window */ 
/* the lower right small window */ 


Loading SMET"); 

*/ 


colorjmenu ■■ defpup("Color %t %F | Red %x0 | Green %xl | Yellow %x2 \ 

1 Blue %x3 I Magenta %x4 | Cyan %x5 | White %x6 | Gray %x7", smet_sat_color); 


orbit_menu « defpup("Define Orbit %t %F | Low Earth | Molniya \ 
I Geosynchronous | Input Parameters", add_satellite); 


add_menu — defpup("New Satellite %t I Define Orbit %m | Define Launch %f \ 
I Read File %f", orbit menu, add launch, ask and read sat file) ; 


param_menu - defpup("Parameter Options %t %F | Classical | Secondary \ 

I Position I Options | Combined", siRetjparam_option); 

picture_menu » defpup("Picture Options %t %F | Color RGB Image (Regular) %xl \ 

I Color RGB Image (Inverted) %x2 | Post Script - Entire Screen %x3 I Post Script \ 

- Main Window %x4", take_picture); 


top_inenu defpup ("SMET menu %t | Add Satellite %m | Maneuver Satellite %f \ 
I Delete Satellite %f | Save Satellites %f | Take Picture %m | Satellite \ 
Parameters %m | Instructions %f | Exit %f ", addjmenu, add_maneuver, 
delete_satellite, ask_and_write_sat_file, picture_menu, param_menu, 
draw instructions, program exit); 
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/* queue the necessary devices */ 

qdevice(MIDDLEMOUSE); 
qdevice(RIGHTMOUSE); 
qdevice(DIAL2); 
qdevice(DIAL3); 
qdevice(DIAL4); 
qdevice(DIALS); 
qdevice(DIAL6); 
qdevice(DIAL?); 

tie(MIDDLEMOUSE,MOUSEX,MOUSEY); 

qreset (); /* initialize the queue */ 

/* dial limits: one minute to 6 hours */ 
setvaluator (DIALS, (short) (oloc)c_speed), 1, 360); 


/* miscellaneous initializations 


*/ 


/* read data files for drawing globe and map */ 

read_world_file(); 

/* Initialize the satellite list */ 

header_sat.next = NULL; 

/* read any sat files listed by user on startup */ 

for (c ■ 1; c < argc; C++) read_3at_file(argv{cj); 

/* Initialize font for windows */ 

fminit () ; 

fmsetfont(fmscalefont(fmfindfont("Times-Roman"), LETTER)); 

/* User Initialition Routine */ 

user init(); 


/* put the different views into the default windows */ 

omni_view_gid — mainwin.gid; 

sat_view_gid — auxwinl.gid; 

merc_view_gid — au::win2 . gid; 

window order - OMNISATMERC; 


/* draw instructions for the user */ 

draw instructions(); 
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/* INIT-WINDOW */ 
/* Set up a window using the data from the given */ 
/* window structure */ 


init_window(win) 

window_ 3 truct *win; /* window structure */ 

( 

/* set the pref position according to structure */ 

prefposition(win->orgx, win->orgx + win—>sizex, 
win->orgy, win->orgy + win->si*ey); 


win->gid — winopen(""); 

/* 

create a window */ 


noborder (); 

/* 

no edges or window manager pop ups 

*/ 

RGBmode(); 

/* 

RGB color specifications 

*/ 

doublebuffer(); 

/* 

no flickering 

*/ 

zbuffer(TRUE); 

/* 

hidden surface removal 

*/ 

lsetdepth(0, 0x7fffff); 

/* 

get the largest Zbuffer resolution 

*/ 

gconfigO ; 

/* 

attach settings to the current window */ 

mmode(MVIEWING); 

/* 

put winod in viewing mode 

*/ 


/* go black and clear both buffers */ 

c3f(Color[BLACK]); 

clear () ; 

swapbuffers () ; 

clear () ; 
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!* ******************************** */ 
/* */ 
/* 3n\et_Jcepler. c */ 

/* by Carlos I. Noriega */ 

/* */ 
/* This file contains the routines to update orbital */ 

/* parameters. It is used by the SMET demonstrator. */ 

/* */ 
/* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/* */ 
/* *********************************/ 


linclude "smet.h" 


/* UPDATE-STATIC-PARAMETERS 
/* See Chapter II of thesis. 

update_static_parameters(sat) 
satellite *sat; 

{ 

/* update the semi-minor axis, period, perigee, apogee when */ 

/* ecc or major have changed due to a maneuver or user input */ 
sat->minor — 3at->major * sqrt(1 - 3at->ecc * 3at->ecc); 
sat->period — pow(sat->major,1.5) * DEG360 * TU; 

3at->apogee_range - 3at->major * (1 + sat->ecc); 

3at->perigee_range - 3at->major * (1 - sat->ecc); 
update_perturb_constants(sat); 

) 


*/ 

*/ 


/* UPDATE-PERTURB_CONSTANTS */ 

/* See Chapter II of thesis. */ 

update_perturb_con3tants(sat) 
satellite *sat; 

{ 

float perturb_constant — J2 / (pow(sat->major, 3.5) * pow((l - sat->ecc * sat->ecc), 2)) 
/* constant for regression of line of nodes */ 

sat->asc_change - - 1.5 * cos(3at->incl) * perturb_constant /TU; 

/* constant for rotation of line of apsides */ 

3at->peri_change - (3 - 3.75 * sin(sat->incl)*3in(sat->incl)) * perturb constant /TU; 
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/* UPDATE-CHANGING-PARAMETERS 
/* See Chapter II of thesis. 

update_changing_parameters(sat, time) 
satellite *sat; 
float time; 

( 

float E, M; /* Eccentric and Mean anomalies */ 

if (clock_running && perturbations) 

( 

/* regression of line of nodes */ 

sat->asc » sat->asc + sat->asc_change * cloclc_speed; 

/* rotation of line of apsides */ 

3at->peri — sat->peri + sat->peri_change * clocJc_speed; 
update_sat_matri:: (sat) ; 

) 

/* use M as initial approximation of eccentric anomaly */ 

sat->E = M = fmod(((time - sat->epoch) * DEG360)/3at->period, DEG360); 

if (M < 0.0) sat->E - M - DEG360 + M; /* ie. time < epoch */ 

/* compute the Eccentric anomaly using Newtonian Approximation */ 

/* normally converges after 3 to 5 iterations at this accuracy */ 
do 
{ 

E = sat->E; 

3at->E « E + (M - E + sat->ecc * sin(E))/(l - sat->ecc * cos(E)); 

) while (nonzerop(3at->E - E)); /* repeat till almost zero */ 

/* compute the true anomaly from the eccentric anomaly */ 

sat->anom acos ( (cos lsat->E) - sat->ecc) / (1 - sat->ecc * cos',sat->E) ) ) ; 

if (sat->E > DEG180) 3at->anom - DEG360 - sat->anom; /* get correct quadrant */ 

/* compute the distance to the satellite */ 

sat->range ■ sat->major * (1 - sat->ecc * cos(sat->E)); 

/* compute the velocity of the satellite */ 

sat->velocity ” sqrt((2/sat->range) - (l/sat->major)); 

/* compute the planar coordinates of the satellite *! 

sat->p — sat->range * cos(sat->anom); 
sat->q ■ sat->range * sin(sat->anom); 

/* compute the XYZ coordinates of the satellite. Don't */ 

/* forget that this is for an IJK system which is */ 

/* different than the display systems XYZ system */ 

3at->y ■ sat->q*sat->matrix[0] (0)+sat->p*sat->matrix[2] (0]+sat->matri::[3] [0); 
sat->z « sat->q*sat->matrix(0][1]+sat->p*sat->matrixl2](l)+sat->matrix(3)(1); 
sat->x ■ sat->q*sat->matrix(0)[2]+3at->p*sat->matrix[2][2]+sat->matrix[3](2); 

/* compute the sat's lat and long for plotting adjustinf for earth spin */ 
sat->lat “ asin(sat->z / sat->range); 

sat->lon — atan2(sat->y, sat->x) - fmod(time * ROTATION_RATE,DEG360); 
if (sat->lon < -DEG180) sat->lon «• sat->lon + DEG360; 

I* save the plot if it is being plotted */ 
if (sat->plot) add_plot(sat); 


*/ 


) 
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/* UPDATE-SAT_MATRIX */ 
/* compute the orbital plane transformation matrix */ 
/* lets let the graphics hardware do some of our */ 
/* tough calculations. what a neat trick! */ 


update_aat_matrix(sat) 
satellite *sat; 

( 


loadmatrix(ident_matrix); 
rot(sat->asc*TO_DEG, 'y'); 

/* 

longitude of the ascending 

node 

*/ 

rot(sat->incl*TO_DEG,'z'); 

/* 

inclination 

of the orbital 

plane 

*/ 

rot(sat->peri*TO_DEG, 'y'); 
getmatrix(sat->matrix) ; 

/* 

argument of 

periapsis 


*/ 


) 


/* ALT_UFDATE-SAT_MATR1X */ 
/* Use for multiprocessing when encorporated. */ 
/* Gets same matri:: as update_sat_matrix. */ 
/* Slightly different than matrix in Chapter 2 of */ 
/* thesis due to difference in coordinate systems. */ 


alt_update_3at_matri:: (sat) 
satellite *sat; 

{ 

float cos^asc = cos(sat->asc), 
sin__asc « sin (sat->a3c) , 
co3__inGl - cos (sat->incl) , 
sin^incl = sin(sat->incl), 
cos__peri = cos (3at->peri), 
3in__peri - sin(3at->peri) ; 


sat->matrix(0)(0] 
sat->matrix[0][1] 
sat->matrix(0](2) 
sat->matrixI1)(0) 
sat->matrix(l][IJ 
3at->matrix(1](2] 
sat->matrix[2][0) 
sat->matrix(2][1] 
3at->matrix[2][2] 
sat->matrix[0)[3] 
sat->matrix[3][0] 
3at->matrix[3)[3] 


-sin_asc*sin_peri + cos_asc*cos_peri*cos_incl; 
cos_peri*3in_incl; 

-cos_asc*sin_peti - sin_asc*cos_peri*cos_incl; 

-cos_a3c*sin_incl; 

cos_incl; 

sin_asc*sin_incl; 

sin_a3c*cos_peri + co3_asc*3in_peri*cos_incl; 
sin_peri*3in_incl; 

cos_a3c*cos_peri - sin_a3C*sin_peri*co3_incl; 
sat->matrix(l] [3] - sat->matrix [2][3] ■ 0.0; 
3at->matrix[3] [1] - sat->matrix[3] [2] - 0.0; 
1 . 0 ; 
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/**yt*******************************/ 
/* */ 

/■* smet_launch. c */ 

/* by Carlos I. Noriega */ 

/* */ 

/* This file contains the routines used by SMET to *! 

/* launch a satellite. */ 

/* */ 

/* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/* */ 

/**********************************/ 

iinclude "smet.h" 

/* APD_LAUMCH */ 

/* Use a series of dialogs to get required info. */ 

/* See Chapter 3 for more information. */ 


add_launch() 

( 

launch_3truct *new_launch - (launch_struct *)malloc(sizeof(launch_struct)) 
char str[15]; 

dialog("Enter name of new satellite.", new_launch->name); 

new_launch->time - ask_time("time at burnout"); 

dialog("Enter the altitude at burnout", str); 

new_launc)i->alt =• a^of (str)/distance_factor [unit_type] ; 

dialog("Enter the velocity at burnout.", str); 

new_launch->v - atof(str)/speed_factor[unit^type]; 

dialog("Enter the latitude at burnout.", str); 

new_launch->lat " atof(str)*TO_RAD; 

dialog("Enter the longitude at burnout.", str); 

new_launch->lon •• atof (str) *TO_RAD; 

dialog("Enter the angle above the horizon.", str); 

new_launch->elev - atof(str)*TO_RAD; 

dialog("Enter the heading angle.", str); 

new_launch->azim - atof(str)*TO_BAD; 

new_launch->next — pending_launches; 

pending_launches — new_launch; 


/* LAUNCH_PENDING_SATS */ 

/* Check all pending to see if it is time to launch. */ 

launch_pending_sats() 

( 

launch_struct "launch, "prev; 

for (launch - pending_launches; launch !— NULL; launch - launch->next) 

{ 

if (launch->time <- clock) 

( 

initialize_orbital_parametera(launch) ; 

if (launch - pending_launches) pending_launches — launch->next; 
else prev->next — launch->next; 
free((char ") launch); 

) 

prev * launch; 

) 
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/* INITIALIZE_ORBITAL_PARAMETERS */ 

/* See Chapter 3 for more information. */ 

initialize_orbital_parameters(launch) 
launch_3truct ^launch; 

( 

char 3tr[70]; 
float p, 

v_radial — launch->v * 3in(launch->elev), 

v_n3 — launch->v * cos(launch->elev) * cos(launch->azim), 
v_ew — launch->v * cos(launch->elev) * sin(launch->azim) 

+ ROTATION_RATE * TU * COS (launch->iat), 

v_hor — sqrt(v_ns*v_ns + v_ew*v_ew); 
satellite *3at; 

make_new_current_3at(launch->name); 
sat — current_sat; 

/* adjust velocity, elevation angle and azimuth for earth spin */ 
launch->v — sqrt (v_radiai*v_radial + v_hor*v__hor) ; 
launch->elev - atan2(v_radial,v_hor); 
launch->azim * atan2(v_ew,v_ns); 

/* compute orbital parameters as discussed in Chapter 3 of thesis */ 

3at->incl “ acos(sin(launch->azim) * cos(launch->lat)); 
sat->range - 1 + launch->alt; 

sat->major - sat->range/(2 - 3at->range * launch->v * launch->v); 
p - 3at->range * v_hor * 3at->range * v_hor; 
sat->ecc - fsqrtd.O - p/sat->major) ; 

sat->anom * atan2(p * tan(launch->elev), p - sat->range); 

3at->E - acos((sat->ecc + cos(sat->anom))/(1 + sat->ecc * cos(sat->anom))); 
if (3at->anom > DEG180) sat->E - DEG360 - sat->E; /* get correct quadrant */ 

sat->epoch » launch->time - pow(sat->major,1.5) * (sar.->E - sat->ecc*sin(3at->E)) * TU 
sat->peti — fmod(atan2(sin(launch->lat),cos(launch->lat) * cos(launch->azim)) 

- sat->anom + DEG360, DEG360); 

3at->asc — fmod(launch->lon + clock*ROTATION_RATE - atan2(sin(launch->lat) 

*sin(launch->azim), cos (launch“>azim) ) + DEG360, DEG360); 

update_static_parameter3(sat); 
update_3at_matrix(sat); 
reset_dial3(); 

if (0.0 > sat->ecc I| sat->ecc >— 1.0) 

( 

sprintf(str,"%s was launched into an open orbit and deleted from the system.", 
sat->name); 
delete satellite(); 

) 

else sprintf(str,"%3 was launched.", sat->name); 
message(str); 
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****^******************i%******** "k / 
/* */ 
/* smet_main.c */ 
/* by Carlos I. Noriega */ 
/* */ 
/* This file main routines used by SMET. */ 
/* */ 
/* Written - September 3, 1990 */ 
/* Modified - September 3, 1990 */ 
/* */ 
/**********************************/ 


finclude "smet.h" 


/* MAIN */ 

/* Initialize everything and enter a loop to check the */ 

/* queue, act on it and update the windowa. */ 

main(argc,argv) 
int argc; 
char **argv; 

K 

initialize_3met(argc,argv); /* windows, menus, queue, variables */ 

while (TRUE) /* loop until we exit via a pop up menu */ 

{ 

handle_queue(); 

user_pre_update_satellites0; /* User defined routine */ 

update_aatellite3(); 

user_post_update_satellites0; !* User defined routine */ 
draw_windows(); 

if (clock_running) clock - clock + clock speed; /• increment clock if on */ 

) 
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/************** This is where I want to do the multi processing ***********/ 


/* UPDATE_SATBLLITES */ 
/* Calling routine for updating all the satellites. */ 
/* Also checks for earth impact and sets Automatic */ 
/* Viewing position if auto_viewer is on. */ 


update_satellites() 

I 

satellite *sat, *previou3_sat; 
float max_apogee “ 1.2; 

previous_sat “ fiheader_sat; 
launch_jpending_sats () ; 

for(sat •• header_sat.next; sat !“ HULL; sat ” sat->next) 

( 

if (sat->maneuver !- NULL £& clock >— 3at->maneuver->time) maneuver_sat(sat); 
u3er_pre_orbital_parameters(sat, clock); 
update_changing_parameters(sat, clock); 
user_post_orbital_parameters(sat, clock); 

if (sat->range <- 1.0) crash_satellite(sat, previou3_sat); 
previous sat “ sat; 

max_apogee = max(3at->apogee_range, max_apogee); 

) 

if (auto viewer) /* put our viewing position out towards Ares */ 

{ 

viewer_x «= max_apogee; 
viewer_y - 1.0; 
viewet_z - 1.0; 

viewer_range ■ sqrt(max_apogee*max_apogee + 2); 

) 

) 


/* DRAW_WINDOWS */ 

/* Calling routine for all the drawing routines. */ 


draw_windows() 
( 


draw_omni_view () ; 

/* 

draw_sat_view(); 

/* 

draw_merc_view(); 

/* 

draw_parameters(); 

/* 

draw_3at_list(); 

/* 


) 


draw the omniscient point of view */ 
draw the satellite point of view */ 
draw the mercator map with plot */ 
show the current satellite's stats */ 
draw the list of satellites */ 
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/* 

it * 

*********************** 

********/ 

!* 



*/ 

/* 


smet maneuver.c 

*/ 

/* 


by Carlos I. Noriega 

*/ 

/* 



*/ 

/" 


This file contains the routines used by SMET 

to */ 

/* 


maneuver a satellite. 

*/ 

/* 



*/ 

/* 


Written - September 3, 1990 

*/ 

/* 


Modified - September 3, 1990 

*/ 

/* 



*/ 

/* 

* * 

*********************** 


♦include "smet.h" 


/* 

ADD 

JMANEUVER 

*/ 

/* 

Use 

a series of dialogs to get required info. 

*/ 

/* 

See 

Chapter 4 for more information. 

*/ 

add 

maneuver () 



{ 

char str(25]; 

man_struct *walker, *maneuver; 

if (current_sat !“ ITOLL) 

( 

maneuver = (man_struct *)malloc(sizeof(man_struct)); 
if (current_3at->maneuver ” NULL) 

( 

dialog("Enter the True Anomoly at the time of the maneuver.", atr); 
maneuver->anom * atof(str)*TO_RAD; 
current_3at->maneuver ■ maneuver; 
compute_maneuver_time(current_sat); 

) 

elae 

{ 

dialog("Enter the number of degrees to travel after the preceding maneuver.", str) 
maneuver->anom — atof(str)*T0_RAD; 

for (walker current_3at->maneuver; walker->next !— NULL; walker — walker->next) 
walker->next — maneuver; 

) 

dialog("Enter (1) for Delta-V/Yaw/Pitch or (2) Velocity/Flight Path/Heading Changes. 

, str) ; 

maneuver->normal_maneuver - atoi(str) % 2; 
if (maneuver->normal_maneuver) 

( 

dialog("Enter the Delta-V for the maneuver.", str); 
maneuver->deltav - atof(str)/3peed_factor(unit_type]; 
dialog("Enter the Pitch for the maneuver.", str); 
maneuver->pitch — atof(str)*TO_RAD; 
dialog("Enter the Yaw for the maneuver.", str); 
maneuver->yaw - atof(str)"TO_RAD; 

I 

else 

( 

dialog("Enter the Velocity Magnitude Change for the maneuver.", str); 
maneuver->deltav - atof(str)/speed_factor[unit_type]; 
dialog("Enter the Flight Path Change for the maneuver.", str); 
maneuver->pitch — atof(str)"TO_RAD; 

dialog("Enter the Heading Change for the maneuver.", str); 
maneuver->yaw - atof(str)*T0_RAD; 

» 90 

maneuver->next - NULL; 

) 
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/* MANEUVER_SAT */ 

/* See Chapter 4 for more information. */ 

maneuver_3at(sat) 
satellite *sat; 

« 

float delta_f_path, delta_heading, new_velocity/ 
old_ecc “ sat->ecc, 
oldjtnajor - sat->major, 
old_incl “ sat->incl, 
old_peri •• sat->peri, 
old_a;ic ■> sat->asc; 

update_changing_parameters(sat, sat->maneuver->time); 

if (3at->maneuver->normal_maneuver) 

( 

compute_deltas(sat, 4delta_f_path, fidelta_heading, 4new_velocity); 
flight_path_change(sat, delta_f_path); 
heading_change(sat, delta_heading); 
velocity_change(sat, new_velocity); 

) 

else 

{ 

flight_path_change(sat, sat->maneuver->pitch); 
heading_change(sat, sat->maneuver->yaw); 

velocity_change(sat, sat->velocity + sat~>maneuver->deltav); 

) 

if (0.0 <- sat->ecc sat->ecc < 1.0) 

{ 

sat->epoch - sat->maneuver->time 

- pow (sat->ma jo*;, 1.5) * (sat->E - sat->ecc * sin(3at->E)) * TU; 
if (sat “ current_sat) reset_dials(); 
update_static_parameters(sat); 
update_sat_matri:: (sat) ; 

) 

else 

( 

message("That maneuver put the satellite into an open orbit and will be ignored."); 

sat->ecc — old_ecc; 

sat->major — old_major; 

sat->incl ” old_incl; 

sat->peri — old_peri; 

sat->asc — old_asc; 

) 

delete maneuver(sat) ; 
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/* COMPUTE_PELTAS */ 

/* See Chapter 4 for more information. */ 

C(m(pute_delta9(sat, delta_f_path, delta_heading, new_velocity) 
satellite *9at; 

float *delta_f_path, *delta_heading, *new_velocity; 

( 

float 

delvp — sat->maneuver->deltav * cos(sat->maneuver->yaw), 
delvps - delvp * sin(aat->maneuver->pitch), 
delvop •• sat->maneuver->deltav * sin(aat->maneuver->yaw), 
vlmodl “ sat->velocity + delvp * cos(aat~>iBaneuver->pitch), 
vlmod2 ■ fsqrt(delvps*delvps + vlmodl^vlmodl); 

*delta_f_path - fatan2(delvps,vlmodl); 

*delta_heading « fatan2(delvop, vlmod2); 

*new_velocity — fsqrt(delvop*delvop + vlffiod2*vlmod2); 

) 

/* COMPUTE_MANEUVER_TIME */ 

/* Set the time of maneuver based on the current time, */ 

/* current anomaly and desired anomaly. */ 

compute_maneuver_time(sat) 
satellite *sat; 

{ 

float c_anom — cos(sat->maneuver->anom), 

E « acos((3at->ecc + c_anom)/(l + sat->eoc * c_anom)); /* Eccentric Anomoly */ 

/* number of revolutions completed plus number of more whole revolutions */ 
int revs - (int)(clock - sat->epoch)/sat->period 

+ (int)(sat->maneuver->anom/DEG360); 

/* go to next time around if we're past it on this revolution*/ 
if (sat->maneuver->anom <- sat->anom) revs++; 

/* make anomoly positive */ 

while (sat->maneuver->anom < 0) sat->maneuver->anom ■ sat->maneuver->anom + DEG360; 

/* correct quadrant */ 

if (fmod(sat->maneuver->anom,DEG360) > OEG180) E » DEG360 - E; 

sat->maneuver->time " sat->epoch + sat->period*(revs + (E - sat->ecc * sin(E))/DEG360); 

) 

/* DELETE_MANEUVER */ 

/* Remove the maneuver from the linked list. */ 

delete_maneuver(sat) 
satellite *3at; 

( 

man_9truct *temp; 

temp ■ 8at->maneuver; 
sat->maneuver - sat->maneuver->next; 
free((char *) temp); 

/* set execution time of next maneuver if there is one */ 
if (sat->maneuver !« MULL) 

( 

sat->maneuver->anom 9at->maneuver->anom + 8at->anom; 
compute_maneuver_time(sat); 92 

) 
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/* FLIGHT_PATH_CHANGE */ 

/* See Chapter 4 for information. */ 

flight_path_change(sat, delta) /* Zn-Plane Maneuver */ 
satellite *sat; 
float delta; 

float p, 

old_anom, 
tan_f_path; 
char strI70]; 

if (nonzerop(delta)) 

I 

tan_f_path ■ sat->ecc * sin(sat->E)/faqrt( 1 - 8at->ecc * 8at->ecc); 
sat->ecc - fsqrt(l - (1 - sat->ecc*sat->ecc> * 

pow((cos(delta) - tan_f_path*3ln(delta)), 2)) ; 

old_anom ■■ sat->anom; 

tan_f_path - (tan_f_path + tan(delta))/( 1 - tan_f_path * tan(delta)); 
p - sat->major • (1 - sat->ecc * sat->ecc); 

sat->anom — fmod(fatan2(p * tan_f_path,p - sat->range) + DEG360, DEG360); 

sat->E « acos((sat->ecc + cos(3at->anom))/(1 + sat->eco * cos(sat->anom))); 

if (sat->anom > DEG180) sat->E — OEG360 - sat->E; /• get correct quadrant •/ 

sat->peri «= fmod (sat->peri + old_anom - sat->anotn + DEG360, DEG360) ; 
sprintf(str, "%3 executed a flight path change of %4.3f degrees.", 
sat->name, delta * TO_DEG); 

message(str) ; 

) 

) 
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/* HEADING_CHANGE *! 

/* See Chapter 4 for information. */ 

heading_change(aat, delta) I* Out-of-Plane Maneuver */ 
satellite *sat; 
float delta; 

( 

float incl_old “ sat->incl, 

u_old - sat->peri + sat->anom, 

u_new — 0.0, 

temp; 

char str(70]; 

if (nonzerop(delta)) 

( 

sat->incl - facos(cos(incl_old)*cos(delta) + sin(inel_old)*cos(u_old)*sin(delta)); 
if (nonzerop(aat->incl)) 

( 

temp * sin(incl_old)*cos(u_old)*cos(delta) - cos(incl_old)*sin(delta); 
u_new — facos(temp/sin(sat->incl)); 

if (sat->y < 0.0) u_new DEG360 - u_new; /* get correct quadrant */ 

) 

sat->peri * fmod(u_new - 3at->anom + DEG360, DEG360); 

3at->asc “ fmod(sat->a3C + fatan2(sin(u_old)*cos(incl_old),cos(u_old)) 

- fatan2(sin(u_new)*cos(sat->incl), cos(u_neH)) + DEG360, DEG360); 
sprintf(str,"%s executed a heading change of %4.3f degreessat->name,delta*TO_DEG); 
message(str); 

> 

) 
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/* VEI,C)CITY_CHANGE */ 

/* See Chapter 4 for information. */ 

velocity_change(sat, new_vel) 
satellite *sat; 
float ne»»_vel; 

{ 

float tan_f_path, cos_fjpath, p, old_anom; 
char str(70]; 

if (nonzerop(new_vel - aat->velocity)) 

( 

if (sat->ecc ” 0) 

if (new_vel > sat->velocity) 

( 

sat->peri — fmod(sat->peri + sat->anom, DEG360); 
sat->E » sat->anom — 0.0; 

) 

else 

( 

sat->peri = fmod(sat->peri + 3at->anom + DEG180, DEG360); 

3at->E — sat->anom « DEG180; 

) 

sat->major - sat->range/(2 - sat->range * new_vel * new_vel); 

tan_f_path « sat->ecc * fsin(sat->E)/sqrt(1 - sat->ecc * sat->ecc); 
co3_f_path — cos (atan (tan_f_r'ath) ) ; 

p » sat->range * 3at->range * new_vel * new_vel * cos_f_path * cos_f_path; 
sat->ecc - f3qrt(1.0 - p/sat->major); 

old_anom ■ sat->anom; 

sat->anom •* fmod(fatan2(p * tan_f_path,p - sat->range) + DEG360, DEG360); 

sat->E “ acos((sat->ecc + cos(sat->arAn))/(1 + sat->ecc * cos(3at->anom))); 

if ( 3 at->anom > DEG180) sat->E - DEG360 - sat->E; /* get correct quadrant */ 

sat->peri = fmod(3at->peri + old_anom - sat->anom + DEG360, DEG360); 

sprintf(str,"%s executed a velocity magnitude change of %4.3f %s.sat->name, 

(nevr_vel - 3at->velocity)*speed_factor[unit_type],speed_unit[unit_type]) 
message(str); 

) 

) 
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/* ******************************** */ 
/* */ 

/* 3met_inap.c */ 

/* by Carlos I. Noriega */ 

/* */ 

/* This file contains the routines to draw the globe */ 

/* used by the SMET demonstrator. */ 

/* */ 

/* Written - September 3, 1990 *f 

/* Modified - September 3, 1990 */ 

/* */ 

/**********************************/ 

finclude "smet.h" 

/* DRAW-MERCATOR *! 

/* */ 


draw_mercator() 

( 

c3f(Color(EARTH]); 

rectf(-DEG180,-DEG90,DEG180,DEG90); /* map background */ 

draw_coa3ts2d(); 
draw_borders2d() ; 
draw_latlon32d() ; 

I 


/* DRAW-COASTS2D */ 

/* Draw the desired data FINE or COURSE. *! 

draw_coa3t32d() 

( 

int c; 

linewidth(2); 
c3f(Color[COAST]); 
if (fineflag) 

for (c - 0; c < FINEPOINTS; C++) 
if (fine2d_data.pen(c]) 

draw(fine2d_data.x[c], fine2d_data.y[c], fine2d_data.z[c]); 
else 

move(fine2d_data.X[c], fine2d_data.y(c], fine2d_data.z[c]); 

else 

for (c - 0; c < COARSEPOINTS; C++) 
if (coarse2d_data.pen(c]) 

draw(coarse2d_data.X[c], coarse2d_data.y(c], coarse2d_data.z[c]) 

else 

move(coarse2d_data.X[cj, coarse2d_data.y[cj, coarse2d_data.z[c]) 

» 


96 



lt_Mp. C 


Hon S«p 3 17:58:59 1990 


2 


/* DRAW-BORDERS2D */ 

/* Draw the border data if desired. */ 

draw_borders2d() 

{ 

int c; 

linewidth(1); 
c3f(Color(BORDER]); 
if (borderflag) 

for (c - 0; c < BORDERPOIKTS; C++) 
if (border2d_data.pen[c]) 

draw(border2d_data.x[c], border2d_data.y[c], border2d_data.z[c]); 
else 

move(border2d_data.x[c], border2d_data.y[c], border2d_data.z[c]); 


/* DRAW-LATLONS2D */ 

/* Draw the latitude and longitude lines on the map. */ 

draw_latlons2d 0 
( 

int c; 

linewidth(1); 

c3f (Color(LATLONG]) ; 

for (c - -6; c <•■ 6; C++) /* draw the longitudes, every 30 degrees */ 

{ 

move(DEG30 * c,-DEG90, 0.0); 
draw(DEG30 * C, DEG90, 0.0); 

) 

for (c “ -3; c <- 3; C++) /* draw the latitudes, every 30 degrees */ 

[ 

move(-DEG180, DEG30 * c, 0.0); 
draw( DEG180, DEG30 * c, 0.0); 

) 
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/* 


****** A 

* * * * * 

^ it it it t 

!* 




*/ 

/* 

smet orbit.c 



*/ 

/* 

by Carlos I. Noriega 



*/ 

/* 




•/ 

/* 

This file contains the routine used by SMET 

to 

*/ 

/* 

draw the orbit ellipses in a viewing window. 

*/ 

/* 




*/ 

/* 

Written - September 3, 1990 



*/ 

/* 

Modified - September 3, 1990 



V 

/* 




*/ 

/* 

itlf'kltititlfititiiititltit'ttliii 

******* 

* * * * * 

•It 

•« 

iinclude "smet.h" 

' 



/* 

DRAW-ORBIT 


*/ 


/* 

draw an orbital ellipse based on 

the global 

*/ 


/* 

parameters: major, minor and 

eccentricity. 

*/ 



draw_orbit(sat) 
satellite *3at; 

{ 

pushmatrix() ; 

multmatrix(sat->matrix); /* set the orbital transformation plane matrix */ 

linewidth(3); /* set the linewidth */ 

c4f(Color[sat->color]); /* go to satellite's color */ 

if (sat->vectors) 

( /* draw the normal vector */ 
move(0.0, 1.0, 0.0); 
draw (0.0, 'S.O, 0.0); 

/* draw the perigee vector */ 

move(0.0, 0.0, 1.0); 

draw(0.0, 0.0, sat->perigee_range); 

) 

/* shift the focus to the origin */ 

translate(0.0, 0.0, -sat->ecc*sat->major); 

/* make an ellipse out of a circle */ 
scale(3at->minor, 1.0, sat->major) ; 

/* put it on equatorial plane */ 

rotate(900,' x'); 

if (sat->fill) 

{ /* draw a simple transparent circle */ 

zbuffer(FALSE); 

blendfunction(BF_SA, BF_MSA); 

circf(0.0,0.0,1.0); 

blendfunction (Br_ONE, BF_ZERO); 

zbuffer(TRUE); 

) 

else circ(0.0, 0.0,1.0) ; /* draw an unfilled oval */ 

popmatrix() ; 

) 
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^'******•***************************0/ 


/* */ 

/* 3inet_parain. c */ 

/* by Carlos I. Noriega */ 

/* * / 

/* This file contains the routines used by SMET to */ 

/* handle parameter changes. */ 

/* */ 

/* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/* */ 


/«***********«***«*************«***/ 
♦include "smet.h" 


smet_param_option(n) /* used by popup menu */ 

i nt n ; 

( 

param_option = n; 

) 


3met_sat_color(n) /* used by popup menu */ 

int n; 

( 

current_3at->color - SAT_RED + n; 

) 


/* 

CHANGE_PARAMETER 



*/ 

/* 

Compute the parameter 

number base 

on the mouse 

*/ 

/* 

position and call the 

appropriate 

change routine. 



change_parameter(y) 
int y; 

int paramnum = (float) (satlistwin.sizey - y) / (float)LINE - 4; 

if (0 <— paramnum && paramnum < 25) 
switch (param_option) ( 

case CLASSICAL: change_classical_param(paramnum) ; brea)c; 

case SECONDARY: change_secondary_param(paramnum); break; 

case POSITIONS: change_position_param(paramnum) ; brea)c; 

case OPTIONS: change_options_param(paramnum); break; 

case COMBINED: change_combined_param(paramnum); break; 

) 

else change_sysparam(paramnum - 29); 

) 
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/* CHANGE_3txxxxxxx_PARAM */ 

/* Change the appropriate parameter. */ 

change_classical_param (paramnum) 
int paramnum; 

( 

switch(paramnum) 

( 

case 0: ask_name(); brealc; 

case 1: as)c__major () ; update_statlc_paraffleter8 (current_8at) ; break; 
case 2: ask_minor_£or_ecc(); update_static_parameter8(current_sat); break; 
case 3: ask_ecc(); update_static_parameters(current_sat); break; 
case 4: ask_incl(); update_sat_raatrix<current_sat); 

updatejperturb_con8tants(current_sat); break; 
case 5: ask_asc(); update_sat_matrix(current_sat); break; 
case 6: ask_peri(); update_aat_niatrix(current_sat) ; break; 

case 7: bad_choice_message(classicalparamblock.params[paramnxmi].name); break; 
case 8: ask_epoch(); delete_plots(current_sat); break; 

) 

if (1 <• paramnum && paramnum <— 6) 

{ 

reset_dials(); 
delete_plots(current_3at); 

} 

) 


change_ 3 econdary_param(paramnum) 
int paramnum; 

{ 

if (paramnum •«— 0) ask_name(); 

else if (paramnum ” 3) ask_epoch(); 
else if (paramnum <•• 6) 

bad_choice_message(secondaryparamblock.paramsIparamnumJ.name); 

) 


change_position__param (paramnum) 
int paramnum; 

if (paramnum 0) ask_name(); 
else if (paramnum <■ 7) 

bad choice_message(positionparamblock.params[paramnum].name); 

) 


change_options_jparam (paramnum) 
int paramnum; 

{ 

char str[50]; 

switch(paramnum) 

( 

case 0: ask_name(); break; 
case 1: dopup(color_menu); break; 

case 2: current_sat->fill -• !current_sat->fill; break; 

case 3: current_3at->vectors “ !current_sat->vectors; break; 

case 4: current_sat->plot ■ !current_8at->plot; deletejplots(current_8at); 
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change_coinbined_param (paramniim) 
int paramnum; 

{ 

i£ (paramnxun <■• 24) switch (paramniun) 

{ 

case 0: ask_naine () ; break; 

case 1: ask_niajor(); update_static_parameters(current_sat); break; 
case 2: ask_minor_for_ecc(); update_static_paraineters(current_sat); break; 
case 3: ask_ecc(); update_static_parameters(current_sat); break; 
case 4: a9k_incl (); update_sat_matrix(current_sat); 

update_perturb_constant3(current_sat); break; 
case 5: ask_asc(); update_sat_matrix(current_sat); break; 
case 6: ask_peri(); update_sat_inatrix(current_sat); break; 
case 8: ask_epoch<); delete_plots(current_sat); break; 
case 21: dopup(color_menu); break; 

case 22: current_sat->fill - !current_sat->fill; break; 
case 23: current_sat->vectors — !current_sat->vectors; break; 

case 24: current_3at->plot ■■ ! current_sat->plot; delete_plots(current_sat) ; break; 
default: bad_choice_ines3age(combinedparamblock.params(paramnum] .name); break; 

] 

if (1 <— paramnum && paramnum <— 6) 

( 

re3et_dials(); 
delete_plot3(current_3at); 

] 

) 
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change_sy3parain (paramnum) 
int paramnum; 

( 

char 3tr{15J; 
satellite *sat; 

switch(paramnum) 

{ 

case 0: clock_running — !clock_running; break; 
case 1: 

clock — ask_time("current time"); 

for (sat - current_3at; sat !- NULL; sat - sat->next) delete_plots (sat); 
break; 
case 2: 

dialog("Enter a clock speed (in minutes per update)",str); 
clock_speed « atof(str); 

setvaluator(DIAL3, (short)(clock_speed), 0, max((short)(clock_speed)"1.2, 360)) 
break; 

case 3: fineflag — ffineflag; break; 
case 4: borderflag = Iborderflag; break; 
case 5: perturbations = (perturbations; break; 
case 6: auto_view€r - !auto_viewer; break; 
case 7; 

dialog("Enter new viewer X position.",str); 
viewer_x = atof(str)/di3tance_factor(unit_type]; 

viewer_range - sqrt(viewer_x"viewer_^x + viewer_y"viewer_y + viewer_z*vieHer_ 2 ); 
auto_viewer -• EALSE; 
break; 
case 8; 

dialog("Enter new viewer Y position.",str); 
viewer_y » atof(str)/distance_factot[unit_typej; 

viewer_range * sqrt(viewer_x*viewer_x + viewer_y*viewer_y + viewer_z*viewer_z); 
auto_viewer « FALSE; 
break; 
case 9: 

dialog("Enter new viewer Z position.",str); 
viewer_z — atof(str)/distance_factor(unit_type]; 

viewer_range — sqrt(viewer_x*viewer__x + viewer_y*viewer_y + viewer_z*viewer_z); 

auto_viewer — FALSE; 

break; 

case 10: unit_type — (unit_type +1) % 3; 

) 

) 


*/ 
*/ 

bad_choice_mes3age(name) 
char "name; 

{ 

char str [50]; 

sprintf(str,"%scannot be changed by the username); 
message(str); 


/* BAD_CHOICE_MESSAGE 
/* 
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/* ASK_xxxxxx - Custom dialog routines used 
/* to ensure the correct input is made, 
float aslc_time (name) 
char *name; 

( 

char msg[75], 

response[20]; 
int hour “ 0, 
min » 0; 
float sec ~ 0.0, 
time; 

sprintf(msg, "Enter the %s.",name); 
dialog(msg, response); 
sscanf(response,"%d%"c%d%"c%f",&hour, &min, &sec); 
if (min >— 0 && min < 60 && sec >~ 0.0 && sec < 60.0) 

if (hour >■ 0 response[0] !“ '-') time « hour"60 + min + sec/60; 
else time - hour*60 - min - sec/60; 
else 
( 

sprintf(msg, "The %s must be entered in 00:00:00 format.", name); 

message(msg); 

time — aa)c_time (name) ; 

) 

return (time) ; 


*/ 

*/ 


ask name() 

( 

char str(45]; 

sprintf (str, "Enter new name for satellite %s current_sat->name); 

dialog(str,current_3at->name); 


a3k_major() 

( 

char str[15]; 

dialog("Enter a Semi-Major Axis.",str); 

current_sat->major - atof(str)/distance_factor[unit_type]; 
if (current_3at->major <— 0.0) 

{ 

message ("Seini-Major Axis must greater than 0."); 
askjma jor () ; 

) 


ask_minor_for_ecc() 

( 

char str[15]; 

dialog("Enter a Semi-Minor Axis.",str); 

current_sat->minor — atof(str)/distance_factor[unit_type]; 

current_sat->ecc - sqrt( 1 - pow(current_sat->minor / current_sat->major, 2)) 
if (current_sat->minor <-0.0 || current_sat->minor > current_sat->major) 

( 

message("Semi-Minor Axis must be between 0 and the Semi-Major Axis."); 
ask minor for ecc(); 

) ■ " ■ 103 
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aak_ecc <) 

( 

char str[15]; 

dialog("Enter an Eccentricitystr); 
current_3at->ecc * atof(3tr); 

if (current_sat->ecc <0.0 |I current_sat->ecc >* 1.0) 

I 

message("Eccentricity must be between 0.0 and 1.0."); 
ask_ecc 0 ; 

) 


ask_incl() 

{ 

char str[15]; 

dialog("Enter an Inclination.",str); 
current_sat->incl — atof(str) * TO_RAD; 

if (current_3at->incl <0.0 || current_sat->incl > DEG180) 

( 

message("Inclination must be between 0 and 180 degrees."); 
as)t_incl () ; 

) 

) 


a3)i_a3C () 

{ 

char str(151; 

dialog("Enter an Ascending Node.",str); 
current_sat->asc - atof(str) * TO_RAD; 

if (current_sat->asc < 0.0 II current_sat->asc > DEG360) 

( 

message("Ascending Node must be between 0 and 360 degrees."); 
as)c_asc () ; 

) 

I 


ask_peri() 

I 

char str[15]; 

dialog ("Enter an Argiiment of Periapsis.", str); 
current_sat->peri - atof(str) * TO_RAD; 

if (current_sat->peri <0.0 || current_sat->peri > DEG360) 

( 

message("Argument of Periapsis must be between 0 and 360 degrees.") 
ask_peri(); 

) 

I 

ask_epoch 0 

I 

curtent_sat->epoch ■ ask_time("epoch"); 

) 
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/**********************************/ 
/* */ 

/ * STnet_queue . c * / 

/* by Carlos I. Noriega */ 

/* */ 

/* This file contains the routines used by SMET to */ 

/* check the queue and handle events not covered */ 

/* in other files. */ 

/* */ 

/* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/* */ 

/********************************★*/ 

finclude "smet.h" 

/* GET_DIAL_VALUE */ 

/* Given a dial and an initial reading from the queue keep */ 

/* reading it until the user stops tuning it. */ 

float get_dial_value(dev,initial_value) 

Device dev; /* which dial we are working with */ 

short initial_value; /* value when first detected on the queue */ 

( 

/* set the current value of the dial */ 

short value — initial value; 


/* read the queue until something other than the */ 
/* same dial is on top or the queue is empty */ 
while (qtest() dev) qread(£value); 

if (dev !- D1AL3) delete_plota(current_sat); 

/* return the last value in float format */ 

return' Afloat)value); 


/* RESET_DIALS */ 

/* Used when a value of a classical parameter is changed */ 

/* without the use of the dials. Keeps everyone happy. */ 

reset_dials() 

I 

/* Set the dials allowing for degree of accuracy desired */ 
setvaluator(DIAL2, (short) (current_sat->peri*TO_DEG*10), 0, 360*10); 
setvaluator(DIAL4, (short)(current_sat->incl*TO_DEG*10), 0, 180*10); 
setvaluator(DIALS, (short)(current_sat->a3c*TO_DEG*10), 0, 360*10); 
setvaluator(DZAL6, (short)(current_sat->major*100), 100, 1000); 
setvaluator(DIAL7, (short)(current_sat->ecc*1000), 0, 999); 
qresetO; /* suppress "setvaluator" induced events. */ 
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I* HANDLE_QUEUE */ 

/* Clears the queue of all pending events. *l 


handle_queue() 
{ 


short 

value; 

/* 

value 

from a device on the queue 

*/ 

while 

(qtest 0) 

/* 

clear 

queue before updating windows 

*/ 


switch(qread(fivalue)) 

( 

case MIDDLEMOUSE: 

if (value 1 ) handlejatLddlemouse (); 
break; 

case RIGHTMOOSE: 

if (value 1 ) dopup (top_menu); 
break; 

case DIAL2; /* Adjust "peri* */ 

if (current_sat !“ NULL) 

{ 

current_sat->peri — get_dial_value(DIAL2,value)*TO_RAD/10; 
update_sat_matrix(current_aat); 

) 

break; 

case D1AL3: /* Adjust "clock speed" */ 

clock_speed - get_dial_value(DIALS,value); 
break; 

case DIAL4: /* Adjust "incl" */ 

if (cutrent_sat !- NULL) 

current_sat->incl - get_dial_value(DIAL4,value)"TO_RAD/10; 
update_sat_matrix(current_sat); 
update_perturb_constants(current_sat); 

) 

break; 

case DIALS: /* Adjust "asc" */ 

if (current_sat !— NULL) 

( 

current_aat->asc " get_dial_value(DIALS,value)"TO_RAD/10; 
update_sat_matrix(current_sat); 

) 

break; 

case DIAL6: /* Adjust "major" */ 

if (current_sat !“ NULL) 

{ 

current_sat->major - get_dial_value(DIAL6,value)/lOO; 
update_static_parameters(current_sat); 

) 

break; 

case DIAL7: /* Adjust "ecc" */ 

if (current_sat !■ NULL) 

I 

current_sat->ecc >■ get_dial_value(DIAL7,value)/lOOO; 
update_static_parameters(currcnt_aat); 

) 

break; 

default: break; 

) /* end of switch statement *l 
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/* PROGRAM-EXIT */ 

/* Exit the program cleanly */ 

program_exlt() 

( 

char str(10]; 

dialog("Would you like to save any satellites before quitting?",str); 

if (str[0J ” 'y' || str[0] — 'Y') ask_and_write_sat_file(); 

message("Thank you for using SMET."); 

winclose(satlistwin.gid); 

winclose(paramwin.gid); 

winclose(auxwinl.gid); 

winclose(auxwin2.gid); 

winclose(mainwin.gid); 

winclose(dialogwin.gid); 

winclose(backwin.gid) ; 

gexitO; /" Close graphics ports */ 
exitO; /* Exit orbital_demo */ 

) 


/* HANDLE-MIDDLEMOUSE */ 

/* Call the appropriate routine based on the */ 

/* mouse position. */ 

handle_middlemou3e() 

( 

short x,y; /* mouse position */ 

qread(4x); /" cet MOUSEX */ 

qread(Sy); /* get MOUSEY "/ 


/* Are we in the area containing the aux windows? */ 

if ( X >— auxwinl.orgx && x <— auxwin2.orgx-fauxwin2.sizex) 

( 

if ( y >“ auxwinl.orgy Sfi y <— auxwinl.orgytauxwinl.sizey ) 
switch_windows(x); 

) 

else if ( X <- paramwin.orgx + paramwin.sizex) change_parameter(y); 
else if ( X >— satlistwin.orgx) switch_current_sat(y); 

) 
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*/ 

*! 

*/ 
*•/ 
*/ 
*! 
*/ 
*/ 
*/ 
*/ 

switch_»r±ndOKS {x) 
short x; 

( 

/* Are we in the left aux window? */ 

if ( X <- auxwinl.orgx -t- auxwinl.sizex ) 

switch (window_order)( 

case OMNISATMERC: 

sat_view_gid ■ mainwin.gid; 
omni_view_gid — auxwinl.gid; 
window_order SATOMNIMERC; 
break ; 

case OMNIMERCSAT: 

merc_view_gid - mainwin.gid; 
onmi_view_gid * auxwinl.gid; 
window_order * MERCOMNISAT; 
break; 

case SATMERCOMNI: 

merc_view_gid - mainwin.gid; 
aat_view_gid ■■ auxwinl.gid; 
window_order “ MERCSATOMNI; 
break; 

case SATOMNIMERC: 

omni_view_gid - mainwin.gid; 
sat_view_gid auxwinl.gid; 
window_order « OMNISATMERC; 
break; 

case MERCOMNISAT: 

omni_view_gid >■ mainwin.gid; 
merc_view_gid — auxwinl.gid; 
window_order » OMNIMERCSAT; 
break; 

case MERCSATOMNI: 

sat_view_gid ■ mainwin.gid; 
merc_view_gid ■ a\ixwinl.gid; 
window_order ■ SATMERCOMNI; 
break; 


/* SWITCH-WINDOWS 

/* See if we are in one of the auxiliaty windows. If we are, 
/* then swap the gid values for the views associated with 

/* the current main window and the auxiliary window which 

/* was selected by the user by pressing the MIDDLEMOUSE. 

/* The window_order variable designates the which views 

/* are in which windows. For exan^le, OMNISATMERC 

/* means that the main window contains the Clnniscient 

/* point of view, auxiliary window 1 contains the 

/* Satellite point of view, and auxiliary window 2 

/* contains the meracator projection. 
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I* Are we in the right aiax window? 
else if (X >- auxwin2.orgx ) 

switch (window_order){ 

case OMMISATMERC: 

mere view_gid “ mainwin.gid; 
oroni~view_gid - a\ixwin2.gid; 
window_order “ MERCSATOMNI; 
break; 

case OMNIMERCSAT: 

sat view_gid >■ mainwin.gid; 
omni_view_gid >■ auxwin2.gid; 
window_order ” SATMERCOMNI; 
break; 

case SATMERCOMNI: 

Omni view_gid " mainwin.gid; 

sat view_gid - auxwin2.gid; 
window order “ OMNIMERCSAT; 
break; 

case SATOMNIMERC: 

mere view_gid ” mainwin.gid; 
sat view_gid - auxwin2.gid; 
window_order ” 1>1ERCC*1NISAT; 
break; 

case MERCOMNISAT: 

sat view gid ” mainwin.gid; 
mere view_gid - auxwin2.gid; 
window_order ” SATOMNIMERC; 
break; 

case MERCSATOMNI: 

orani view_gid = mainwin.gid; 
merc_view_gid *= auxwin2.gid; 
window_order OMNISATMERC; 
break; 

) 
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! * ******************************** */ 


/* */ 

/* smet_sat.c */ 

/* by Carlos 1. Noriega */ 

/* */ 

/* This file contains the routines used by SMET to */ 

/* draw the satellite in a viewing window. It was borrowed */ 

/* from the ASAT program written by Jim Zanoli and Do Kyeing */ 

/* Ok. It was modified and converted to procedures. *! 

/* */ 

/* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/* */ 


!**********************************/ 
finclude "smet.h" 


/* DRAW-SATELLITE */ 
/* draws a satellite in the correct attitute in */ 
/* the appropriate spot in the orbital plane */ 
/* as set by the global orbital parameters */ 


/* anom, r_p and r_q. The satellite is scaled */ 

I* with respect to the major axis of the orbit.*/ 

draw_satellite(sat) 
satellite *sat; 

( /* scale size increases with distance */ 

float sat_acale - 0.005 * viewer_range; /* maybe vr*vr ?? */ 

short angle; /* dummy loop counter for making the ball */ 

pushmatrix(); 

multmatrix(sat->matrix); /* set the orbital transformation plane matrix */ 

linewidth(2); /* set the linewidth */ 


translate(sat->q,0.0,sat->p); 
rot(sat->anom*TO_DEG, 'y'); 
rotate(900,'X'); /* 

scale(sat scale,sat scale,sat 


/* move it in the orbital plane */ 

/* point towards the earth 
lay it on its side */ 

scale); /* make much smaller than the earth 


c3f(Color[WHITE]); /* 

go 

white 

/* draw left wing */ 
pushmatrix 0 ; 
scale(3.0,0.7,0.1) ; 

/* 

make 

translate(-1.2,0.0,0.0); 

/* 

move 

draw_cube (); 
popmatrix (); 



/* draw right wing */ 
pushmatrix 0 ; 
scale(3.0,0.7,0.1); 

/* 

make 

translate(1.2,0.0,0.0); 

/* 

move 

draw_cube(); 



popmatrix(); 



c3f(Color[sat->color)) ; 


/* 


*/ 


it long, flat, and tl...n */ 
it out to the side */ 


it long, flat, and thin */ 
it out to the side */ 


make the body the right color */ 
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/* draw satellite body */ 
draw_cylinder(); 

/* draw satellite middle globe part */ 
pushmatrix 0 ; 

translate(0.0,-3.5,0.0); /* move below the main cylinder */ 

c3f(Color[YELLOW]); /* go yellow */ 

/* draw a 0.8 radius ball */ 

for(angle—0;angle<“3 600;angle^angle+2 0 0) 

( 

rotate(300, 'x'}; 
circf(0.0,0.0,0.8); 
rotate(300, 'y') ; 
circf(0.0,0.0,0.8) ; 
rotate(300,' z') ; 
circf(0.0,0.0,0.8); 

) 

popmatrix () ; 

c3f(Color[3at->color]) ; 

/* draw front small body 
pushmatrix() ; 
translate(0.0,-5.0,0.0) ; 

3cale(1.0,0.3,1.0) ; 
draw_cylinder0 ; 
popmatrix() ; 

popmatrix() ; 


/* make the body the right color 

*/ 

/* move below the ball */ 

/* make it shorter */ 


*/ 
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/* DRAW-CYLINDER */ 

/* draw a 6 unit long by 1 unit radius cylinder */ 

/* centered at the origin and sitting upright */ 

draw_cylinder() 

{ 

short angle; 

/* draw the sides of the cylinder */ 

/* by drawing 72 filled rectangles */ 
pushmatrix(); 

for (angle**0; angle<—3600; angle—angle+50) 

( 

rotate(50,'y'); 
rectfi(-1,-3,1,3); 

) 

popmatrix (); 

/* draw the top of the cylinder */ 

pushmatrix() ; 

translate(0.0,3.0,0.0); /* move it to the top */ 

rotate(900,'x'); /* lay it flat */ 

circfi (0,0,1); 
popmatrix(); 

/* draw the bottom of the cylinder */ 
pushmatrix(); 

translate(0.0,-3.0,0.0); /* move it to the bottom */ 

rotate(-900,'x'); /* lay it flat */ 

circfi (0,0,1); 
popmatrix(); 
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/* DRAW-CUBE */ 

/* draw a cube, used for reshaping to any box shape */ 

draw_cube() 

( 

/* draw the front */ 
pushmatrlx(); 
translate(0.0,0.0,1.0); 
rectfi(-1,-1,1,1); 
popmatrix(); 

I* draw the back */ 
pushmatrlx(); 
translate(1.0,0.0,0.0); 
rotate(900,'y'); 
rectfi(-1,-1,1,1); 
popmatrix 0 ; 

/* draw the right */ 
pushmatrlx(); 
translate(0.0, 0.0, -1.0) ; 
rotate(1800,'y') ; 
rectfi(-1, -1,1,1) ; 
popmatrix() ; 

/* draw the left */ 
pushmatrlx(); 
translate(-1.0,0.0,0.0) ; 
rotate(-900,' y' ) ; 
rectfi(-1,-1,1,1) ; 
popmatrix() ; 

/* draw the top */ 
pushmatrlx(); 
translate(0.0,1.0, 0.0) ; 
rotate(-900,' x' ) ; 
rectfi (-1,-1,1,1) ; 
popmatrix(); 

/* draw the bottom */ 
pushmatrlx(); 
translate(0.0,-1.0,0.0); 
rotate(900,'x'); 
rectfi(-1,-1,1,1); 
popmatrix(); 
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/* */ 

/* 9met_u3er.c */ 

/* by Carlos I. Noriega */ 

/* */ 

/* This file contains the user routines called by SMET. */ 

/* They currently do nothing. */ 

/* */ 

/* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/* */ 

/******************«***************/ 


finclude "smet.h" 

U3er_init() 

\ 

/* ex, read sat files, init variables etc. */ 

\ 


U9er_new_3at() 

I 

/* ex. modify new current_3at */ 

) 


U3er_pre_update_3atellite9() 

( 

/* ex. check for exit, compare sat positions, maneuver sats, launch sats, etc. */ 

) 

user_post_update_satellite3 <) 

{ 

/* ex. check ..v^r exit, compare sat positions, maneuver sats, launch sats, etc. *! 

» 


u9er_pre_orbital_parameter3(sat, time) 
satellite *3at; 
float time; 

( 


/* e.x. perturb orbit etc. */ 

) 

user post_orbital_parameters(sat, time) 
satellite *9at; 
float time; 

{ 

/* ex. perturb orbit etc. */ 

» 
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/* ******************************** */ 
/* */ 
/ * 3inet_windows. c * / 

/* by Carlos I. Noriega */ 

/* */ 
/* This file contains the routines used by SMET to */ 

/* draw all the windows. */ 

/* */ 
/* Written - September 3, 1990 */ 

/* Modified - September 3, 1990 */ 

/* */ 
/********************************* *^ 


finclude "smet.h" 


/* DRAW-OMNI-VIEW */ 

/* Draw the globe, orbital ellipse, satellite and */ 

/* and the UK axis at the center of the window */ 

/* form an omniscient point of view. */ 

draw_omni_view() 

( 

satellite *3at; 

/* select the proper window and and clear it */ 

winset(omni view_gid); 

c3f(Color[BLACK]); 

zclear 0 ; 

clear (); 

/* make the box large enough to hold any possible rotation */ 

/* of the orbital plane. Keep the globe in the center. */ 

ortho(-viewer_range, viewer_range, -viewer_range, viewer_range, 0, 100); 

loadmatrix(ident_matrix); 

/* draw earth background */ 

pushmatrix () ; 

look=*t(0.0, 0.0, viewer_range, 0.0, 0.0, 0.0, 0.0); 
c3f(Color[EARTH]); 
circf(0.0,0.0,1.0) ; 
popmatrix 0 ; 

/* place our viewer */ 

lookat(viewer_y, viewer_z, viewer_x, 0.0, 0.0, 0.0, 0.0); 

/* draw the orbital ellipses */ 

for(sat ■■ current_3at; sat !“ NULL; sat — 3at->next) draw_orbit(sat); 

/* draw the UK a::es */ 

draw_axes(); 

/* draw the earth w/ rotation */ 

draw_globe (ROTATION_RATE * clock); 

/* draw the satellites */ 

for(sat - current_3at; sat !“ NULL; sat -> sat->next) draw_3atellite(sat) 
swapbuffers0; /* show it */ 
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/* DRAW-SAT-VIEW */ 
/* Show that portion of the globe that ia visible *! 
/* to the satellite. This is best represented */ 


/* by a "fish-eye" view given by using perspective. */ 

draw_sat_view{) 

\ 

float del_z; 

/* select the proper window and clear it */ 

winset(sat_view_gid); 

c3f(Color[BLACK]) ; 

zclear(); 

clear (); 

/* if we're above the surface draw it */ 
if (current_sat !- NULL £& current_aat->range > 1.001) 

( 

/* compute the fov based on the range and make the "/ 

/" box deep enough to handle the near (viewable) */ 

/" half of the globe. */ 

del_z — l/current_sat->range; 

perspective(asin(del_z)"3790/PI, 1.0, (current_sat->range -1)/2,current_3at->range) 
loadmatrix(ident_matrix); 

pushmatrix(); 

lookat(0.0, 0.0, current_sat->range, 0.0, 0.0, 0.0, 0.0); 
translate(0.0, 0.0, del_z); 

/* draw a ring around the globe the color of the current satellite */ 

c3f(Color{current_3at->color]); 

circf(0.0, 0.0, 1.02 " 3qrt(l - del_z*del_z)); 

/* draw earth background */ 
c3f(Color [EARTH]) ; 

circf(0.0, 0.0, sqrtd - del_z"del_z)) ; 
popmatrix() ; 


/* look from the satellite to the center of the earth */ 

/* remember that UK system not same as IRIS' XYZ system */ 

lookat (current_3at->y, current_sat->z, current_sat->x, 0.0, 0.0, 0.0, 0.0); 

/* draw the earth w/ provided rotation */ 
draw_globe(ROTATION_RATE " clock); 

) 

swapbuffers0; /* show it */ 
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/* DRAW-MERC-VIEW */ 
/* Draw a mercator projection map and plot all */ 
/* lat/long 


draw_merc_view() 

I 

satellite *3at; 

plot_struct ^plotter; 

int k; /* dummy loop counter */ 

/* get the right window */ 

winset(merc_view_gid); 

/• clear out any viewing or projection matrices left */ 
/* on the window and set up a simple 2D screen */ 

ortho2(-DBG180-0.001,DEG180+0.001,-DEG120,DEG120); 
loadmatrix(ident_matrix); 

c3f(Color[BLACK]); /* go to black and clear window */ 

clear(); 
zclear () ; 


draw_mercator(); 

/* 

draw 

the map 

*/ 

draw_3un(); 


draw 

golden circle for sun 

*/ 


for(sat = current_3at; sat != NULL; sat ” 3at->next) if (3at->plot) 

( 

c3f(Color[3at->color]); 

for(plotter «• 3at->plots; plotter !■ NOLL; plotter - plotter->next) 
dtaw_dot(plotter->lon, plotter->lat); 

» 

swapbuffers0; /* show it */ 
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/* DRAW-DOT */ 

/* Used in drawing the Mercator map. */ 

/* Draw a dot at the given x,y coordinate. Since it */ 

/* is a small dot, it is faster to draw an 8 point */ 

/* circle than the standard circf circle. */ 

draw_dot(x,y) 

float X, y; /* coordinates of center */ 

raove(x, y, 1.0); 
rpmv(-0.00873,-0.02107, 0.0) ; 
rpdr(0.01234,0.01234,0.0) ; 
rpdr(0.0,0.01745,0.0) ; 
rpdr(-0.01234,0.01234,0.0) ; 
rpdr(-0.01745,0.0,0.0) ; 
rpdr(-0.01234,-0.01234,0.0) ; 
rpdr(0.0,-0.01745,0.0) ; 
rpdr(0.01234,-0.01234,0.0) ; 
pclos 0 ; 

) 


/* DRAW-SUN */ 
/* Used in drawing the Mercator map. */ 
/* Draw a transparent ellipse to represent the sun's */ 
/* position. */ 


draw_sun() 

{ 

float sun_lat - SUN_INCL * sin(clock/SUN_TU), 

Sun_lon - fmod (clock * (1/SUN_TU - ROTATION_RATE) - DEG180, DEG360); 

c4f(Color[SUN]); 

zbuffer(FALSE) ; 

blendfunction (BF_SA, Br_MSA); 

circf(sun_lon, sun_lat,1.0); /* western hemisphere */ 

circf(sun_lon + DEG360, sun_lat, 1.0); /* eastern hemisphere */ 
blendfunction(BF_OME, Br_ZERO); 
zbuffer(TRUE); 
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DRAW-PARAMETERS */ 

/* Go through the appropriate parameter block and print */ 

/* each parameter in the right apot in the parauneter window. */ 

draw_j>arameter3 () 

{ 

ParBlock *3at_block; 
int k; 

winaet(paramwin.gid); 
c3f(Color[BLACK]); 
clear (); 

ortho2(0,paramwin.3izex,0,paramwin.3izey); 
loadmatrlx(ident matrix); 


/* current satellite parameter block */ 
I* dumny loop counter */ 

/* go to the parameter window */ 


/* draw satellite parameters */ 

c3f(Color[WHITE]); I* write a title */ 

cmov2i(10, paramwln.sizey - 3 * LINE); 
fmprstr("Satellite Parameters:"); 


if (current_sat !» NULL) 
( 


switch (param_option) { 


case CLASSICAL: 
case SECONDARY; 
case POSITIONS: 
case OPTIONS; 
case COMBINED: 


sat_block 
sat_block 
sat_block 
sat_block 
sat block 


Sclassicalparamblock; break; 
isecondaryparamblock; break; 
Spositionparamblock; break; 
&optionsparamblock; break; 
Scombinedparamblock; break; 


) 

for (k - 5; k < 3at_block->numparam3 +5; k++) 

( 


) 


cmov2i(10, paramwin.sizey - k " LINE); 

draw_parameter(3at_block->param3[k - 5]rcurrent_3at->color, 

(int)current 3at-(int)&header 


) 


sat) 


/" draw system parameters "/ 

c3f(Color[WHITE]); /" write a title "/ 

cmov2i(10, paramwin.sizey - 32 "LINE); 
fmprstr("System Parameters:"); 

for (k • 34; k < sysparamblock.numparams + 34; k++) 

( 

cmov2i(10, paramwin.sizey - k * LINE); 
draw_parameter(sysparamblock-patams[k - 34]rWHITE,0); 

) 

swapbuffers0; /* show it */ 


119 






aaMit_window«. c Mon S«p 3 18:31:43 1990 8 

/* DRAW-PAFAMETER */ 

/* Draw the contents of a parameter record in the */ 

/* desired color. Offset is used to get the value from */ 

/* the desired satellite. */ 

draw_par«imeter (param, color_code, offset) 

Par param; /* a parameter record */ 

int color_code, I* what color to draw the data */ 

offset; /* offset of the string pointer */ 

K 

c3f(Color[WHITE]); 
fmprstr (param. name); 

param.value 7 (char *)((int)param.value -f offset); 
c3f(Color[color_code]); 

switch (param.type) { /* choose the correct output routine */ 

case DISTANCE: draw_distance(param.value); break; 

case SPEED: draw_speed(param.value); break; 

case FLOAT: draw_float(param.value); break; 

case DEGREE: draw_degree(param.value); break; 

case LAT: draw_lat(param.value); break; 

case LON: draw_lon(param.value); break; 

case TIME: draw_time(param.value); break; 

case STRING: fmprstr(param.value); break; 

case COLORCODE: fmprstr(color_names(*(int *)param.value]); break; 
case UNITCODE: fmprstr(unit_namea(*(int *)param.value]); break; 

case FLAG: if(*(int *)param.value) fmprstr("On"); else fmprstr("Off"); break; 

) 


/* DRAW-xxxxxxx "/ 

/♦ Customized formatting routines that coerce the value */ 

/* into the format we want. */ 

draw_distance(val) 
char *val; 

( 

char buf(15]; 


sprintf(buf, "%4.2f %a", 
fmprstr(buf) ; 


"(float *)val*distance_factor[unit_type], 
distance_unit(unit type]); 


draw_speed(val) 
char *val; 

( 

char buf[15]; 


sprintf(buf, "%4.2f %s", 
fmprstr(buf); 


*(float *)val*3peed_factor[unit_type), 
speed_unit[unit_type]); 
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draw_float(val) 
char *val; 

( 

char buf[10]; 

sprintf(buf, "%3.2f", *(float *)val); 
fmprstr(buf); 


draw_de 9 ree(val) 
char *val; 
t 

char buf[10ij; 

sprintf(buf, "%3.2f", *(float *)val*TO_DEG); 
finprstr (buf); 


draw_lat(val) 
char *val; 

( 

char buf[10]; 

if (* (float *)val >■■ 0.0) 

sprintf (buc*, "%3.2f N", fabs(*(float *) val) *TO_DEG); 
else 

sprintf(buf, "%3.2f S", fabs(*(float *)val)*T0_DEG); 
fmprstr(buf); 


draw_lon(val) 
char *val; 

( 

char buf[10]; 

if (* (float *)val >- 0.0) 

sprintf(buf, "%3.2f E", fabs(*(float *)val)*TO_DEG); 
else 

sprintf(buf, "%3.2f W", fabs(*(float *)val)*TO_DEG); 
finprstr (buf) ; 


draw_time(val) 
char *val; 

( 

char buf[10]; 

int hour •» * (float *)val / 60, 

min - * (float *)val - hour * 60, 

sec ■ (*(float *)val - (hour * 60 + min)) * 60; 

if (*(float *)val >- 0 || *(float *)val <- -60) 

sprintf(buf, "%02d:%02d:%02d", hour, abs(min), abs(sec)); 

else sprintf(buf, "-0:%02d:%02d", abs(min), abs(sec)); 

fmprstr(buf); 


) 
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/* DRAW-SAT-IilST */ 

I* List out the linked list of aatellitea */ 

/* vertically in the satlistwin. */ 

draw_aat_li 3 t() 

I 

aatellite *aat; 

int y « aatliatwin.aizey - 3 * LINE; /* firat line of text */ 
wlnaet(aatliatwin.gid); 

c3f(Color[BLACK]); t* clear the window */ 

clear(); 

ortho 2 ( 0 ,aatliatwin.aizex, 0 ,aatliatwin.aizey); 
loadmatrix(ident_matrlx); 

cinov 2 i ( 20 , y) ; 


c3f(Color[WHITE]); 

/* 

write a title 

*/ 

fmpratr("SMET Satellite Liat: 
y - y - LINE; 

") ; 

/* 

blank line 

*/ 

for(3at - current sat; aat !- 
( 

y - y - LINE; 
cmov2i (25,y) ; 
c3f(Color[3at->color]) ; 

NULL; 3at = sat->next) 


/* 

current sat color 

*/ 

fmpr3tr (aat->naine) ; 

» 

awapbuffer3 (); 

/* 

write a name 


/’» 

show it 

*/ 
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/* DRAW_INSTRUCTIONS */ 

/* Loop through the instructions array and draw it on the */ 

I* Main window, then wait for the user to press and */ 

/* release the rightmouse button. */ 

draw_instructions() 

( 

int k; /* dummy loop variable */ 

winset(mainwin.gid); /* get the main window */ 

c3f(Color[BLACK]); 

clear () ; 

zclear (); 

c3f(Color[CYAN]); 

/* clear out any viewing or projection matrices left */ 

/* on the main window and set up a simple 2D screen */ 
ortho2(0,mainwin.sizex,0,mainwin.sizey); 
loadmatrix(ident_matrix); 

for (k *0; k < 28; k++) 

( 

cmov2(75,700 - k * LINE); /* move to the left margin */ 
fmprstr(instructions[k]); /* draw the line of text */ 

) 

swapbuffers0; /* show it */ 

while (qcead(sk) I* RIGHTMOUSE); /* wait for press */ 
while (qread(Sk) != RIGHTMOUSE); /* wait for release */ 
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aak«flltt Mon S«p 3 17:59:01 1990 1 

f/* *********************************/ 
#/* */ 

#/* makefile for SMET */ 

#/* by Carlos 1. Noriega */ 

#/* */ 

#/* Written - September 3, 1990 */ 

#/* Modified - September 3, 1990 */ 

#/* */ 


#/* *********************************/ 

# coR^iler switches 
CC - cc 

INCLUDE - -I. 

CFLAGS - -O -C -D$(LANGUAGE) 

LIBS - -Imalloc -Ifm -lgl_s -Im \ 

imagesupport/libnpsimage.a ~zyda/utah/lib/librle.a -limage 
SRCS — smet_addelate.c\ 

3met_axes.c\ 
smet_dialog. c\ 

3met_flies.c\ 
smet_globe.c\ 

3met_init. 

3met_kepler.c\ 

3met_lc.u: ,h,c\ 

3met main.c\ 

3met_maneuver.c\ 

3met_map.c\ 

3met_orbit.c\ 

3met_param.c\ 
smet_queue.c\ 

3met_3at.c\ 

3met_user.c\ 
smet_windows.c\ 

3met_var3.c 

OBJS — 3met_addelete.o\ 

smet_axe3.o\ 

3met_dialog.o\ 

3met_file3.o\ 

3met_globe.o\ 
smet_init.o\ 
smet_kepler.o\ 
smet_launch.o\ 

3met_main.o\ 

3met_maneuver.o\ 

3met_map.o\ 

3met_orbit.o\ 

3met_parara. o\ 
smet_queue.o\ 
smet_sat.o\ 

3met_user.o\ 
smet_windows.o\ 

3met_var3.o 

3met: $(OBJS) 

cc -O $(INCLUDE) $(OBJS) -o smet $(LIBS) 

clean: 

rm -f $ (ALL) 

Oecho "make clean new con^leted" 

neat: 

rm -f .[BC]* .em* *.o *.BA* *.CK* core 

9echo "make neat new completed" {24 


( 




. c. o: 


$(CC) $(CFLAGS) $(INCLUDE) $< 




LIST OF REFERENCES 
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